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The  trend  in  communication  systems  is  towards  processing  larger  bandwidth  signals.  As 
the  bandwiddi  of  the  signal  increase,  the  fundamental  operations  on  the  signal  become  more 
computationally  intensive.  Optical  signal  processing  systems  have  a  natural  parallel  structure,  so 
that  many  confutations  can  be  done  simultaneously,  therefore  reducing  the  signal  processing  time. 
For  example,  a  1024x1024  point  Fourier  transform  of  a  complex  signal  can  be  computed  optically 
in  the  time  it  takes  light  to  traverse  a  firaction  of  a  meter,  or  about  1  nanosecond. 

As  the  spectral  environment  becomes  more  congested,  wide  bandwidth  systems  are 
frequently  comipted  by  narrowband  interference.  The  narrowband  interference  could  be 
intentional,  or  it  could  be  some  FM  radio  energy  corrupting  a  wideband  satellite  message.  We 
desire  to  excise  the  narrowband  interference  in  order  to  recover  the  intended  signal.  Since  the 
optical  Fourier  transform  exist  as  a  spatial  light  distribution  that  is  physically  accessible,  we  can 
use  spatial  light  modulators  to  physically  block  the  undesired  light,  or  narrowband  interferers,  in 
the  Fourier  transform  plane.  In  optical  signal  processing  we  use  acousto-optic  cells  to  spatially  and 
temporally  modulate  light  in  real-time  to  handle  signals  having  bandwidths  that  range  from  50  MHz 
to  2  GHz. 

In  the  early  1980’s  Erickson  [1-3]  performed  analyses  and  experiments  dealing  with  signal 
excision.  His  main  interest  was  in  finding  the  notch  shape  and  depth  as  various  parameters  of  his 
optical  signal  processor  were  varied.  He  used  wires  in  the  frequency  plane  as  spatial  light 
modulators  to  perform  the  excision.  In  1980  Lee  et  al.  [4]  reported  on  experiments  in  which  they 
used  a  threshold  detector  array  to  control  spatial  light  modulators  in  the  frequency  plane.  They 
used  a  PLZT  ceramic  electro-optic  modulator  in  which  18  dB  of  excision  was  attained,  and  23  dB 
when  an  opaque  light  block  was  used.  In  1982  Roth  [5]  reported  that  he  used  an  array  of 
wideband  photodiodes  in  the  frequency  plane  to  perform  notch  filtering;  a  maximum  notch  depth  of 
35  dB  was  attained.  Anderson  et  al.  [6]  were  the  first  to  demonstrate  excision  using  a  GaAs 
photodetector  array.  They  obtained  a  40  dB  excision  depths,  therefore  making  the  array  useful  for 


correlation  operations  in  optical  ackqjtive  filters. 

Bandstettcr  and  Grieve  [7-8]  built  an  optical  notching  filter  that  uses  a  reoirsive  technique 
to  obtain  die  final  notch  depth.  They  used  a  prograimnable  spatial  filter  in  which  liquid  crystal 
fingers  block  the  undesiied  light  Their  experiments  involves  passing  an  RF  pulse  through  their 
optical  excisor  and  examining  the  effect  their  system  has  on  die  pulse;  they  also  added  narrowband 
inteiferers  to  the  pulse,  excised  them  using  dieir  recursive  architecture,  and  then  examined  the 
output  pulse  at  RF. 

The  Presort  Processor  is  an  optical  signal  processor,  developed  by  the  Harris  Corporation, 
to  excise  narrowband  interference  from  wideband  signals  using  frequency  plane  excision.  The 
Presort  Processor  is  located  at  Rome  Air  Development  Center,  where  it  is  maintained  in  the 
Photonics  Center  The  frequency  plane  excision  is  performed  by  a  pair  of  multi-channci  acousto- 
opdc  devices  currently  controlled  by  a  bank  of  swin:hes.  We  have  designed  and  constructed  a 
detection,  processing,  and  control  system  so  that  die  notching  can  eventually  be  automated. 

The  major  objective  of  our  work  is  to  modify  the  Presort  Processor  so  that  the  frequency 
excision  is  open-loop  adaptive.  Our  major  accomplishments  are: 

•  the  design,  construction,  and  testing  of  a  detection  system  that  obtains  the  Fourier 
transform  of  a  wideband  signal 

•  the  design,  construction,  and  testing  of  a  system  that  adaptively  searches  a  wideband 
spectrum  for  narrowband  interferers. 

•  an  analysis  of  area  detection  in  the  Fourier  plane  versus  point  detection  in  the  image  plane 
of  optical  processors. 

•  die  effect  of  excising  narrowband  interference  from  wideband  signals  is  studied  by 
simulating  the  removal  of  symmetric  frequency  con^nents  of  a  pulse  at  an  intermediate 
frequency,  and  analyzing  the  pulse  shapes  at  baseband;  these  computer  simuladons  are 
confirmed  experimentally  by  utiliang  the  Presort  Processor: 
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•  the  operatioD  of  the  Presort  Processor  is  demonstrated  by  adding  a  narrowband  inteiferer 
to  a  pulse,  excising  the  inteifeier  at  an  intermediate  frequency,  and  examining  the 
distortion  on  the  pulse  at  baseband. 

In  Chapter  2,  we  give  a  discussion  of  optical  signal  processing  fundamentals;  this  includes 
(1)  a  discussion  of  wideband  signals  and  advantageous  processing  architectures,  (2)  a  derivation  of 
the  optical  Fourier  transform,  (3)  a  brief  discussion  of  spatial  light  modulators  with  insight  into 
desirable  real-time  spatial  li^t  modulators,  and  (4)  an  example  of  an  optical  spectrum  analyzer  In 
Chapter  3,  we  derive  the  output  photocurrent  for  die  opdcal  architecture  on  which  the  Presort 
Processor  is  based.  We  then  show  that  area-detection  in  the  frequency  plane  is  equivalent  to  point- 
detection  in  the  image  plane.  In  Chapter  4,  we  describe  the  Presort  Processor  in  detail,  design  a 
system  that  searches  for  narrowband  interferers  in  wideband  spectrums,  and  give  experimental 
results  of  the  systera  In  Chapter  5,  we  present  computer  simulations  and  experimental  results  of 
pulse  experiments  where  we  used  the  Presort  Processor  to  notch  various  frequency  components. 

In  Chapter  6,  we  present  our  conclusions  and  suggestions  for  future  work  in  this  area. 
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2.  Optical  Signal  Processing  Fundamentals 

In  diis  chapter  we  explore  some  of  the  fundamental  concepts  of  optical  signal  processing. 
First,  a  discussion  of  signal  characteristics  is  presented  with  some  insight  into  real-time,  large 
bandwidth  processing  systems.  Second,  we  develop  die  optical  Fourier  transform.  Third,  we 
examine  spatial  light  modulators  and  discuss  which  ones  are  useful  for  real-time  signal  processing. 
Finally,  we  provide  an  example  of  an  optical  spectrum  analyzer 
2.1  Discussion  of  signals 

The  trend  in  communication  systems  is  towards  processing  larger  bandwidth  signals.  As 
the  bandwidth  of  the  signal  increase,  the  fundamental  operations  on  the  signal  become  more 
computationally  intensive.  One  measure  of  the  complexity  of  a  signal  is  its  time  bandwidth 
product  The  time  bandwidth  product  is  given  by  TW,  where  T  and  W  are  the  duration  and 
bandwidth  of  the  signal,  respectively.  The  minimum  number  of  samples  needed  to  accurately 
represent  any  finite  duration,  band  limited  signal  is  Ns2TW,  based  upon  the  Nyquist  sampling 
rate.  Therefore,  a  signal  that  has  a  shon  duration  and  a  narrow  bandwidtii  is  less  conq>lex  than  a 
signal  that  is  of  long  duration  and  large  bandwidth  (i.e.  N  is  larger  for  the  latter  case). 

Computers  have  yet  to  provide  real-time  processing  of  signals  tha:  have  large  time- 
bandwidth  products.  To  process  a  signal  which  has  a  time-bandwidth  product,  TW,  in  real-time, 
TW  operations  must  be  coii9>leted  in  Tq  seconds,  where  To=l/2W  to  satisfy  the  Nyquist  sampling 
rate.  Therefore  TW/ro=2TW2  operations  must  be  con^leted  to  process  signals  in  real-time.  A 
parallel  architecture  is  advantageous  for  processing  con^lex  signals  in  real-time;  the  time  required 
to  conq}lete  a  particular  task  is  significantly  reduced  by  doing  many  operations  in  parallel  Optical 
signal  processing  systems  have  a  natural  parallel  structure,  so  that  many  coirputations  can  be  done 
optically  that  are  cumbersome,  or  impossible,  to  do  digitally.  For  example,  a  1024x1024  point 
Fourier  transform  of  a  complex  signal  can  be  conputed  optically  in  the  time  it  takes  light  to  traverse 
a  firaction  of  a  meter,  which  is  about  1  nanosecond. 
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The  spectrum  of  a  signal  often  contains  important  infoimadon  about  Ae  signal.  The 
Fourier  transform  provides  a  way  of  obtaining  the  spectrum  of  a  signal.  The  continuous  time 
Fouria  transform  integral  is  given  by 

-where  F(f,t)  is  the  frequency  spectrum,  evaluated  at  time  t  over  the  past  T  seconds  of  sig^ 
history,  and  f(z)  is  the  signal.  Notice  that  die  spectrum  is  computed  over  a  sliding  time  interval  of 
duration  T.  We  represent  the  Fourier  transfoim  in  this  form  because  the  spectrum  of  a  signal  needs 
to  be  computed  “on  die  fly”  for  real-time  signal  processing.  In  real-time  systems  the  signal  slides 
through  the  time  window,  T,  so  that  an  instantaneous  spectrum  is  computed.  Figures  2.1a  and 
2.  lb  show  a  representation  of  the  sliding  time  window  and  the  corresponding  instantaneous 
spectrum^ 

We  give  a  simple  example  to  gain  further  insight  into  how  real-time  systems  operate.  In 
Figure  2.1a  we  show  a  signal  that  is  received  from  a  radio  station  in  which  the  message  is  spread 
over  a  wide  bandwidth  using  a  spread  spectrum  technique  before  transmission;  e.g.  the  radio 
message  could  be  a  song  named  “Nothing  Qimpatcs  2  U”  by  Sinead  O’Conner  The  song  is 
distorted  due  to  narrowband  interference  encountered  during  transmission.  To  listen  to  hex  song, 
we  need  to  compute  the  spectrum  of  the  signal,  and  remove,  or  excise,  the  narrowband  interferers. 

For  this  process  to  be  accomplished  in  real-time,  the  spectrum  of  the  signal  must  be 
continuously  displayed  and  the  narrowband  interferers  excised  as  they  appear  Therefore,  the 
spectrum,  shown  in  Figure  2.1b,  of  a  section  of  the  signal  (T  seconds)  is  computed.  The 
narrowband  interferers  are  excised,  using  spatial  light  modulators,  and  the  signal  is  then 
demodulated  to  obtain  the  message.  We  therefore  need  a  system  that  displays  the  spectrum  of  a 
signal  instantaneously,  we  now  show  how  the  optical  Fourier  transform  meets  our  needs  for  real¬ 
time  processing. 

2,2  The  optical  Fourier  transform 

Figure  2.2  shows  a  system  that  computes  the  optical  Fourier  transform.  For  convenience 


we  set  die  lens  one  focal  length  away  from  the  object  and  the  Fourier  planes.  Let  f(x),  which  is  m 
plane  P^,  be  illuminated  by  monochromatic  light  of  wavelength  X ;  note  that  f(x)  could  be  recorded 
on  photographic.film  or  any  one  of  a  number  of  spatial  li^t  modulators.  First  we  calculate  die 
li£^t  distribution  in  plane  P2  due  to  die  signal  f(x).  From  Goodman  [9],  we  know  that  the  light 
distribution  at  plane  P2  is  the  Fresnel  transform  of  f(x).  given  by : 

8(w)-j  f(x)e^^  ^  dx, 

•''i  (2.2) 

where  die  scaling  factors  of  die  Fresnel  transform  have  been  ignored.  We  know  that  a  positive 

spherical  lens  is  represented  by  the  function, 

j-^w^ 

g(w)  =  e  ^  ,  (2.3) 

where  F  is  the  focal  length,  and  X  is  the  wavelengdi  of  the  li^t  source.  Therefore,  die  light 

distribution  at  plane  P3  is  given  by. 


m(w)  ss  g(w)h(w)  s  ^  f  f(x)e  dx. 

(2.4) 

Expanding  the  exponential  and  simplifying  m(w)  leads  to 

~  I  f(x)e  ^  dx. 

(2.5) 

We  find  the  light  distribution  at  plane  P4  by  invoking  the  Fresnel  transform  on  m(w)  from  plane  P3 


(2.4) 


to  P4 : 


<"1 

Substituting  equation  (2.5)  into  (2.6)  we  get, 

r  1 


I  /(x)e  ^  dx. 


U  ‘  J 

Expanding  and  rearranging  equation  (2.7)  we  find  that 

and,  by  combining  terms  in  the  exponentiaL  we  write  equation 


dx. 

(2.6) 

:  e  ^  dw. 

(2.7) 

dxdw^ 

(2.8) 

/ 

dxdw. 


(2.9) 
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We  rewrite  die  exponential  term  to  conqilete  the  square  widi  respect  to  the  w  variable  to  get 


dxdw. 


We  change  the  (»der  of  integration  and  expand  the  exponential  to  reveal  that 


XF 


Pi 


F«)=f./Wl 

Equation  (2.1 1)  is  rearranged  as, 

f«)  =  J  /U) 


dw 


~J-^-x^-^^+x^+4^-2x^) 

^  dx. 


dw 


.lit  . 


(2.10) 


(2.11) 


(2.12) 


This  result  is  result  is  similar  to  the  Fourier  transform,  except  for  the  bracketed  term.  Hence,  to 
put  equation  (2.1 1)  into  standard  form,  we  need  to  show  that  the  bracketed  term  is  not  a  function 
of  ^  or  X.  We  do  not  go  throu^  die  analysis  here;  die  interested  reader  can  find  a  more  corrqilete 
analysis  by  VanderLugt  [10]. 

Hence,  if  the  various  scaling  factors  are  ignored,  the  light  distribution  in  plane  P4 ,  due  to  a 
function  illuminated  by  coherent  light  in  plane  P^,  is  given  by. 


F(.i)  =  l  nx'ie'^^dx. 


(2.13) 


Notice  that  this  is  a  one-dimensional  Fourier  transform  integral,  where  is  the  spatial 
frequency  and  x  is  the  independent  variable.  Furthermore,  we  note  that,  since  x  is  a  measure  of 
length,  the  spatial  frequency  has  units  of  cycles  per  unit  length. 

The  two-dimensional  Fourier  transform  integral  is  a  direa  extension  of  equation  (2.13)  and 
is  given  by, 

J  ^^dxdy. 

(2.14) 

Equation  (2.14)  is  useful  when  dealing  widi  two-dimensional  spatial  light  modulators,  such  as 
two-channel  orthogonal  acousto-optic  cells,  or  with  optical  image  processing. 

2.3  Spatial  light  modulators 

We  have  discussed  some  of  the  basic  ideas  involved  with  optical  signal  processing.  To 
process  signals  optically  we  need  a  device  that  interfaces  electrical  signals  and  optical  signals  (i.e. 
light  distributions).  When  we  need  to  convert  an  analog  signal  to  a  digital  signal,  we  use  an  analog 
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to  digital  converter;  similarly,  when  converting  digital  to  analog  signals,  we  use  a  digital  to  analog 
convertet  In  optics  we  use  a  spatial  light  modulator  to  convert  an  electrical  signal  to  an  optical 
signal. 

Spatial  light  modulators  need  the  following  characteristics  to  process  real-time  signals. 

First,  they  need  a  large  time  bandwidth  product  and  wide  dynamic  range;  the  large  space 
bandwidth  product  ensures  that  complicated  signals  can  be  processed.  As  with  most  systems  in 
signal  processing,  it  is  advantageous  to  have  linear  characteristics  over  a  given  range  so  that  linear 
systems  theory  ^>plies.  A  rapid  response  time  is  also  needed  so  the  spatial  li^t  modulators  can  be 
updated  and  refreshed;  if  die  response  time  is  slow,  real-time  operation  is  not  possible.  Good 
efBciency,  phase  control,  and  geometric  fidelity  are  odier  characteristics  which  are  desirable  for  a 
spatial  light  modulator 

The  performance  of  acousto-optic  spatial  light  modulators  is  extremely  weU  suited  for  use 
in  real-time  operations.  In  1921  BriUouin  [11]  predicted  how  light  and  sound  waves  would 
interact  However,  it  was  not  until  1932  diat  Debyc-Sears  [12]  and  Lucas-Biquard  [13]  obtained 
experimental  results.  Three  years  later,  in  1935,  Raman  and  Nath  [14]  mathematically  modeled  the 
interaction  between  light  and  sound  waves.  It  is  worth  noting  that  these  experimental  results  were 
obtained  without  the  use  of  a  good  coherent  light  source  (e.g.  a  laser),  which  makes  these 
discoveries  even  more  fascinating.  In  the  late  1950’s  and  early  1960’s,  the  laser  was  developed 
and  paved  the  way  for  better  use  of  spatial  light  nxxlulators.  In  the  1970*s,  some  35  years  after 
Raman  and  Nath  had  successfully  modeled  the  phenomenon,  spatial  light  modulators  with  good 
optical  quality  and  large  bandwidths  were  developed. 

2,3.1  Acousto-optic  spatial  light  modulators 

In  this  section  we  give  a  brief  discussion  of  how  acousto-optic  cells  operate.  Consider  the 
^ousto-optic  cell  shown  in  Figure  2.3.  Sound  waves  are  used  to  modulate  the  li^t  in  a  linear 
•fashion  over  a  given  range.  These  sound  waves  are  in  water,  glass,  or,  as  more  commonly  used 
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today,  an  exodc  crystal  A  piezoelectric  transducer  is  mounted  onto  one  end  of  the  material  and  an 


electrical  signal  is  ^lied  to  die  transducer  The  signal  creates  an  acoustic  wave  in  die  material 

inducing  different  pressures  in  the  material.  These  pressure  changes  in  the  material  force  the 

density  to  change,  which  causes  a  change  in  the  index  of  refraction.  The  index  of  refraction 

changes  correspond  to  the  ^plied  voltage.  In  turn,  the  optical  path  lengdi  changes  in  the  direction 

of  that  the  li^t  propagates,  and  the  light  is  phase  modulated. 

Acousto-optic  cells  have  bandwiddis  that  range  from  SO  MHz  to  2  GHz,  depending  on  the 

interaction  material  used.  Since  die  acousto-optic  cells  cannot  support  low  frequencies,  the  signal 

usually  is  mixed  with  a  local  oscillator  signal  so  diat  the  signal  falls  into  the  cells  frequency  range. 

The  length  of  the  cell  with  which  die  light  interacts  is  Z,  and  L  is  the  length  of  the  cell  In 

using  acousto-optic  cells  there  are  two  primary  mottes  of  operation,  the  Bragg  mode  and  the 

Raman-Nath  mode.  To  distinguish  between  the  two  modes  we  define  a  Q  factor 

^  2jrAZ 

'*0  ^  (2.15) 

where  k  is  the  frequency  of  lighi  Oq  is  the  index  of  refraction  of  the  material,  and  A  is  the  acoustic 
wavelength  inside  the  acousto-optic  cell  If  Q>10  the  acousto-optic  cell  is  operating  in  the  Bragg 


mode  and  if  Q<10  the  cell  is  operating  in  the  Raman-Nath  mode.  Note  that 

J  (2.16) 

where  v  is  the  acoustic  velocity  of  the  waves  inside  the  cell  and  f  is  the  frequency  i^lied  at  the 
transducer  We  substitute  equation  (2.16)  into  equation  (2.  IS)  yielding, 

G=— iz4- 

"0  (2.17) 

Therefore,  the  Q  factor  is  proportional  to  the  square  of  the  applied  frequency  of  the  signal.  We 
increase  the  qiplied  frequency  of  the  acousto-optic  cell  to  pass  from  one  mode  of  operation  to 
anodier. 


The  Raman-Nath  mode  of  operation  yields  several  diffracted  orders  due  to  die  lack  of 
interaction  length  of  the  light  with  the  acoustic  waves.  For  most  applications,  we  block  all  but  one 
diffracted  order,  u.sing  only  a  small  portion  of  the  incident  light  as  depicted  in  Figure  2.4. 
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■  Since  good  efficiency  is  desired,  die  Bragg  mode  of  operation  is  generally  used.  Figure  2.5  shows 
'die  Bragg  mode  of  operation  which  yields  the  undifEracted  light  and  one  diffiacted  order  Since, 
efficiencies  for  the  Bragg  mode  can  be  as  hi^  as  90%,  this  mode  is  used  in  practical  real'time 
systems. 

2.3.2  Angle  of  deflection 

Let  us  now  examine  how  the  drive  frequency,  spatial  frequency,  diffiaction  angle,  and 
-  acoustic  wavelength  are  related.  Figure  2.6  shows  an  acousto^optic  cell,  widi  the  various 
parameters  displayed  and  how  they  are  related.  Frcxn  equation  (2. 16)  we  know  how  the  temporal 
frequencies  relate  to  the  acoustic  wavelength  within  the  acousto-optic  cell  The  spatial  fequency  is 


defined  as 


1 

a  =  — . 
A 


(2.18) 


We  substitute  equation  (Z16)  into  equation  (2.18)  to  get  that 

a.L. 

V  (2.19) 

Hence,  for  every  drive  frequency  f  there  is  a  spatial  frequency  corresponding  to  the  tcrr^ral 

frequency  of  the  signal  The  physical  angle  0  is  related  to  a  by  the  relationship  that 

V  (2.20) 

This  result  further  further  illustrates  how  the  light  is  deflected  proportional  to  the  applied  frequency 
of  the  cell,  Le.,  as  the  frequency  of  the  signal  increases  the  angle  of  deflection  increases. 

Therefore,  the  acousto-optic  cell  converts  the  signal  into  an  angular  spectrum. 

2.4  Simple  spectrum  analysis 

Figure  2.7  shows  an  optical  spectrum  analyzer  We  assume  that  the  acousto-optic  cell  is  in 
the  Bragg  mode  of  operation,  and  that  the  lens  is  one  focal  lengtii  away.  A  signal  f(t),  applied  to 
the  transducer  of  the  acousto-optic  cell  becomes  a  traveling  pressure  wave  inside  the  cell  to 
produce  a  function  of  space  and  time,  represented  by 

\  V  2/  (2,21) 

where  T/2  is  a  time  delay  factor  due  to  the  optical  axis  being  in  the  middle  of  the  acousto-optic  cell. 
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We  assume  that  the  signal  is  already  in  the  passband  of  the  celL  The  signal  completely  fills  the  cell, 

whose  length  L  is  the  product  of  the  velocity  of  the  acoustic  wave  and  the  transit  time  of  die  cell  T: 

LsvT.  The  li^tdistribution  after  die  acousto-opdc  cell  is  given  by, 

sU,/)  =  jma{x)f{t  -  T  /  2  -  X  /  v).  (2.22) 

where  a(x)  is  die  aperture  function,  and  m  is  the  modulation  index. 

The  Fourier  transforming  property  of  the  lens  shows  that  at  plane  P2,  we  get 

F(a,t)  =  f^f^i 

where  the  integradon  is  done  on  the  traveling  wave  over  die  endrety  of  the  cell,  and  the  aperture 
funcdon  is  rectangular,  we  have  ignored  the  moduladon  index  to  illustrate  the  basic  idea.  We  make 


a  variable  substitution  by  letting 

X  T 

u^t - —, 

V  2 

(2.24) 

so  that 

dx 

du  =  ~. 

V 

(2.25) 

Hence  F(a,t)  becomes 

f  (a,f)  s 

(2.26) 

We  rearrange  equation  (2.26)  yielding, 

{ 

F(a,t)  =  J fiu)ej^^du. 

t-T 

(2.27) 

Equadon  (2.27)  is  the  Fourier  transform  of  f(u),  where  the  dme  window  is  T  seconds  in  duration. 
As  the  signal  f(t)  con^letely  fills  the  cell,  we  recognize  that  the  most  recent  part  of  f(t)  is  at  the 
transducer,  while  the  end  of  the  acousto-optic  ceil  holds  the  part  of  the  signal  that  occurred  T 
seconds  ago.  This  explains  why  the  integradon  is  equadon  (2.27)  is  over  a  total  time  interval  of  [t- 
(t-T)]=T  seconds. 

The  spectrum  that  ^ipears  in  plane  P2  is  called  the  instantaneous  spectrum  of  f(t)  for  a 
couple  of  reasons.  We  recall  that  the  most  recent  part  of  the  signal  is  at  the  beginning  of  the 
acousto-opdc  cell  and  the  signal  is  continuously  fed  into  the  cell.  The  integradon  is  done  over 
time,  therefore  the  spectrum  of  the  T  seconds  of  die  signal  in  the  cell  appears  instantaneously  in 
plane  P2,  because  the  spectrum  is  computed  in  die  time  it  takes  light  to  travel  less  than  a  meter. 


(2.23) 


which  is  about  one  nanosecond  Since  the  acousto-opdc  cells  hold  a  maximum  of  50 
microseconds  of  die  signal,  the  spectrum  of  the  signal  is  available  instantaneously. 

From  this^liscussion  we  recognize  that  ui  optical  system  can  be  implemented  to  display  the 
instantaneous  spectrum  of  a  signal  The  time  window  of  die  Fourier  transform  integral  slides 
along  the  signal  which  is  a  desirable  quality  because  we  can  find  the  Fourier  transform  of  the  signal 
instantaneously  as  the  signal  is  received  Hence,  if  acousto-opdc  cells  are  used  as  the  spatial  light 
modulators,  the  spectrum  of  wideband  signals  can  be  confuted  in  real-time.  A  more  con^lete 
discussion  of  spectrum  analysis  is  given  in  Chapter  3. 


3.  Signal  Excision 
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In  this  chi^tter  we  discuss  the  Mach-Zehndcr  inteiferometer  when  used  as  an 
interferometric  ^ectium  analyzer  We  derive  the  heterodyne  transfoim  and  give  an  overview  of 
frequency  plane  notch  filtering.  We  compare  inu^  plane  detection  to  Fourier  plane  detection,  and 
show  diat  the  photocunents  for  Fourier  plane  detection  and  image  plane  detection  using  an  area 
photodetector  are  equal.  We  show  that  the  cross-product  photocurrent,  for  image  plane  detection 
using  a  point  photodetector,  is  equal  to  the  cross-product  photocurrent  for  Fourier  plane  detection 
using  an  area  photodetector,  but  that  the  bias  photocunents  are  not  equal  We  analyze  the  effect 
diat  the  two  detection  schemes  have  on  the  signal-to-noise  ratio.  Finally,  we  give  a  discussion  of 
die  effect  that  the  reference  beam  position  has  on  notch  depth. 

In  Section  2.5  we  derived  an  example  of  a  power  spectrum  analyzer  in  which  a  CCD 
photodeteaor  array  is  used  to  detect  the  power  spectrum  of  the  signal  The  photodetector  array  is  a 
square  law  device  so  that  the  phase  of  the  signal  important  in  some  signal  processing  applications, 
is  lost  One  way  of  preserving  die  phase  of  a  signal  in  optical  signal  processing  is  to  use 
heterodyne  detection. 

To  perform  heterodyne  detection,  we  add  a  reference  beam  to  the  signal  beam,  and  collect 
the  light  onto  a  wideband  photodetector.  Figure  3.1  shows  an  interferometric  spectrum  analyzer 
with  a  signal  acousto-opdc  cell  in  the  bottom  branch  of  die  inteiferometer,  and  a  lens,  used  to 
condition  die  reference  beam,  in  the  top  branch.  In  the  following  exanqiles  the  reference  beam  is 
expressed  in  terms  of  its  response  in  plane  P2-  The  signal  and  reference  beams  combine  at  the 
beamcombincr,  and  are  Fourier  transformed  by  lens  L4  to  plane  P4,  where  they  are  collected  onto  a 
wideband  photodetector. 

3.1  The  heterodyne  transform 

Figure  3.1  shows  an  acousto-optic  cell  in  the  lower  branch  of  the  interferometer  that 
receives  the  electrical  signal,  f(t),  and  modulates  the  light  in  both  space  and  time.  The  light  in  the 


14 

upper  branch  of  the  interferometBr  passes  throu^  a  lens  at  plane  to  create  a  point  source  at 

■  plane  P2.  The  point  source  occurs  at  exactly  the  same  distance  from  the  beam  combine  as  does  the 
acousto-optic  ceU.  The  signal  and  reference  beams  combine  at  the  beam  combiner,  and  are  Fourier 


transfarmed  by  lens  L4  onto  a  wideband  ar photodetector  located  in  plane  P4. 

The  baseband  signal  s(t)  is  translated  into  the  bandpass  frequency  response  of  the  acousto- 
opdc  cell  by  using  a  local  oscillator,  cos(2}tfct).  The  positive  diffracted  order,  directly  after  the 
acousto-optic  cell,  is  represented  by 

\  2  vj  (3.1) 

where  m  is  the  modulation  index,  and  a(x)  is  the  £q)erture  function.  We  find  the  Fourier  transform 
of  the  signal  beam  to  be 

F;(a,r)=  r 

J—  y  2  vj  (3.2) 

where  the  Fourier  transform  is  performed  with  respea  to  the  input  spatial  variable  x.  We  remove 

the  factors  that  are  not  a  function  of  x  from  the  integral  to  get 


'+  (a, t)  =  pne  ^ ^ J  ^ j(x)5^r  “  y  “  £ j^>2jor(a-Of 


(3.3) 


We  now  mm  our  attention  to  the  reference  beam,  which  is  a  function  of  x  only.  Since  a 
point  source  is  created  a  distance  of  Xq  from  the  optical  axis  in  plane  P2,  die  reference  beam  is 


represented  by 


r(x)  =  sin 


X-Xq 


(3.4) 


where  do  is  die  size  of  the  spot,  as  measured  from  the  centroid  of  the  sine  function  to  the  first  zero. 
We  Fourier  transform  the  reference  beam  to  find  that 


/?(a)  =  rect 


(X^<Xr  jlxa  X. 
- -  e 


where  aco=l/2do,  and  where  the  exponential  term  is  generated  by  the  off-axis  spot  Also  note  that 
we  are  interested  only  in  the  spatial  region  where  the  spectrum  of  the  signal  and  the  reference  beam 
are  located,  which  is  in  the  vicinity  of  a^. 


In  opucal  processors  the  light  can  be  detected  in  either  the  image  plane  or  the  Fourier  plane, 


15 


depending  on  which  one  is  the  most  convenient  In  Sections  3.1.1  and  3.1.2  we  show  that  the 
ouq>ut  cross-product  photocurrents,  when  the  light  is  detected  in  the  Fourier  plane  and  image 
plane,  lespectivdy,  are  equal. 

3.1.1  Detection  in  the  Fourier  plane 


(3.7) 


We  find  the  light  intensity  of  the  combined  beams: 

/(a,0  =  +l/?(a)p +2Re[F;(a,/)/?*+(a)].  ^^.6) 

A  large  area  photodetector  integrates  the  intensity  ova*  the  output  plane  such  duit  the  photodetector 

current,  g(t),  is  given  by 

=  (3  7) 

where  P(a)  is  the  photodetector  aperture  function.  In  this  analysis  we  assume  that  P(a)  is  uniform 

over  the  surface  of  the  photodetectoc  Since,  the  intensity  bias  terms  firom  equation  (3.6)  arc  not  a 

function  of  time,  they  are  rejected  by  a  bandpass  filter  so  that  the  cross-product  photocurrent  is 

given  by 

giit)  *  2Rc[^j^F+(a,/)/?*(a>iaj. 

We  substitute  equations  (3.3)  and  (3.5)  into  equation  (3.8)  to  get 

*3(0  =  2R. 


We  evaluate  the  term  inside  the  parentheses  by  integrating  over  a  to  give 


’  j2x(a^+a„Xx-x„)  _  J2it(,a^-a„y(,x-x„y 
=  j2K{x-x„)  • 


We  simplify  equation  (3.10)  to  get 

F{a)  = 


n{x-Xg) 


Now  recall  that  d()3l/2aco  to  give 

F{a) 


■*»^sinc 

<  ) 


We  substitute  equation  (3.12)  into  equation  (3.9)  to  get 
ft(t)  =  2Re 


(3.9) 


(3.10) 


(3.11) 


(3.12) 


x-x, 


-  ja(x)s|^f  -  Y  - 


(3.13) 
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Rearranging  equation  (3.13),  we  get  that 
1  . 

—pne  '• 

d„ 


g3(r)*2Rd 


Combining  the  e;q>onential  terms  inside  the  integral  we  get 


(3.14) 


g3(0  =  2Re 


(3.15) 

The  sine  function  is  a  suitable  form  of  a  delta  function  because  no  signal  information  is  altered  by  it 


in  the  frequency  domain.  We  therefore  apply  the  sifting  property  of  the  delta  fiir^on  to  get 


g3(r)  =  2Re 


•aix-Xo)5 


i;-H) 


(3.16) 


We  combine  terms  and  perform  the  teal  operation  to  give, 

where  g3(t)  is  proportional  to  the  input  signal,  s(t),  with  an  associated  time  delay.  We  realize  diat 
g3(t)  is  in  the  same  form  as  the  input  double  sideband  modulation  signal,  but  with  some  scaling 
factors.  We  note  that  a(xo)  is  the  aperture  function  evaluated  at  the  equivalent  probe  position  Xq,  f^. 
is  the  same  frequency  as  the  input  carrier,  and  xq/v  is  a  time  delay  which  is  dependent  on  the 
position  of  the  probe.  Furthermore,  to  recover  the  baseband  signal  s(t)  we  need  to  synchronously 
demodulate  g3  (t).  Since  we  have  shown  that  the  signal  can  be  recovered,  we  therefore  can  broaden 
the  range  of  operation  and  inclement  arbitrary  filter  functions,  by  controlling  the  shape  of  the 
reference  beam  r(x),  oar  inplementing  a  frequency  plane  filter  function  R(a),  both  of  which  can 
be  functions  of  time. 

3.1,2  Detection  in  the  image  piane 

In  this  architecture,  we  detected  the  light  in  the  Fourier  plane  of  the  acousto-optic  cell  to 
obtain  the  signal,  however,  we  now  show  that  we  get  the  same  result  if  die  light  is  detected  in  the 
image  plane.  Figure  3.2  shows  the  interferometric  spectrum  analyzer  with  a  point  photodetector  in 
the  image  plane  of  the  acousto-opdc  cell  and  the  image  plane  of  plane  P2.  Since  the  photodetector 


is  in  the  image  plane  of  the  acousto-optic  cell,  we  represent  the  signal  beam  at  the  photodetector 


reversed  coordinates: 


(3.18) 


where  the  various  parameters  are  die  same  as  before.  Similaiiy,  we  represent  the  reference  beam  at 


the  photodetector  as 


(3.19) 


The  photodetector  integrates  the  light  intensity  over  its  surface,  and  the  bias  terms  are  rejected  by  a 


bandpass  filter  The  cross-product  pbotocurrent  term  is  given  by 

=  2Re^|]‘^4(tt,r)r*+(u)duj. 

We  substitute  equation  (3.18)  and  (3.19)  into  equation  (3.20)  to  get 


(3.20) 


g3(f)  =  2Re 


The  sine  function  is  approximated  as  a  delta  function  to  give 


(3.21) 


g3(f)  =  2Re  ;me  ^  ^a(-a)5^r— je  ^''^5{-u-XQ)du  . 


(3.22) 


We  use  the  sifting  property  of  the  delta  function  and  perform  the  real  operation  to  get 

Except  for  the  reversal  of  die  spatial  coordinates,  equation  (3.23)  is  equal  to  die  cross-product 
photocuirent  represented  by  equation  (3.17).  The  signal  is  therefore  propagating  in  the  opposite 
direction  dian  in  the  acousto-optic  cell  and  which  does  not  affea  the  ouq>ut  detection  process. 
Using  this  architecture,  we  can  therefore  detea  the  light  in  either  the  frequency  plane  or  the  image 
plane  of  the  acousto-opdc  cell. 

3.2  Notch  filtering  in  the  frequency  plane 

Wideband  communications  systems  are  frequently  corrupted  with  narrowband  interference. 


The  narrowband  interference  could  be  intentional  or  it  could  be  the  radio  FM  band  corrupting  a 
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satellite  transmission.  In  any  event,  tfie  interference  corrupts  die  message  at  the  receiver:  We 

desire  to  excise  the  narrowband  interference  to  assist  in  the  recovery  of  the  intended  message. 

Figure  33  shows  a  wideband  spectrum  that  is  corrupted  by  three  narrowband  interferers. 

As  discussed  before,  the  Fourier  transform  exist  as  a  spatial  light  distribution  that  is  directly 

accessible.  We  can  therefore  use  spatial  light  rindulators  to  physically  block,  or  notch  filter,  the 

flpHftdnihle  light  in  the  Fourier  transform  plane.  Figure  3.4  shows  a  notching  filter  function  that 

excises  the  narrowband  interferers.  For  notch  filtering,  the  ideal  spatial  light  modulator  is  one  that 

has  an  infinite  on/off  ratio  and  sufficient  resolution  to  notch  any  width  interference  signal  that  is  at 

jiy  position  relative  to  the  spatial  light  riKxuilator  elements. 

We  now  derive  the  general  form  of  the  cross-product  photocurrent  for  fiequency  plane 

notching.  The  Fourier  transform  of  liic  signal  is  represented  by  F(a,t),  the  Fourier  transform  of 

the  reference  beam  is  R(a),  the  fiequency  plane  filter  is  H(a),  and  the  photodetector  response  is 

P(a).  We  collect  all  Ac  light  onto  a  wideband  photodetector,  filter  the  bias  terms,  and  retain  the 

cross-product  term  to  give 

g3(r)  =  2 Re| £*^F{a,f)R*(a)|//(af  |F(a)|^t/a|. 

Equation  (334)  gives  a  general  form  for  the  cross-product  photocurrent  when  an  arbitrary  filter 

function  H(a)  is  placed  in  the  fiequency  plane. 

As  an  example,  we  let  the  reference  beam  at  the  Fourier  plane  be 

(3.25) 

where  R(a)  is  the  Fourier  transform  of  a  point  source  located  off-axis  by  a  distance  Xq.  If  we 


assume  a  cw  interferer  at  a  fiequency  fj,  the  Fourier  transform  of  the  signal  is 

where  ACa-Oj)  is  the  Fourier  transform  of  the  illumination  function. 
g3(/)  =  2Re' / 


(3.26) 


(3.27) 


If  we  sweep  the  interferer  fiequency  over  the  bandwidth  of  the  system,  we  get  a  convolution  of 
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A(a)=smc(aL)  with  the  magnitude  squared  of  die  filter  function  multiplied  by  the  Fourier 
transform  of  die  reference  beam;  P(a)  defines  the  sweep  range  in  terms  of  spatial 
frequencies.  AsBrickson  [1]  showed,  we  get  some  interesting  results  as  we  let  xq  vary. 

Figure  3.5a  shows  the  notch  filter  and  A(a)  before  the  fiequency  is  swept  over  the  band  and 
Hgure  3.5b  shows  the  response  afterwards.  The  sidelobe  structure  of  the  sine  function  shown  in 
Figure  3.5a  is  caused  by  the  illumination  of  the  signal  acousto-optic  cell  and,  mote  fundamoitally, 
the  finite  length  of  the  cell  To  excise  all  of  the  energy  at  fj,  we  need  to  block  the  mainlobe  plus  all 
of  the  sidelobes.  Blocking  all  the  sidelobe  energy  is  not  practical  because  that  would  remove  too 
much  signal  energy.  We  therefore  note  that  the  notch  depth  is  a  function  of  the  apoture  function, 
which  imposes  some  fundamental  limits  on  notch  depth. 

3.3  Problem  description  of  point  detection  versus  area  detection 

We  now  explore  the  differences  between  point  detection  in  the  image  plane  and  area 
detection  in  the  Fourier  plane.  Figure  3. 1  shows  an  interferometric  spectrum  analyzer  with  a  lens 
in  the  upper  branch  and  an  acousto-optic  cell  in  the  lower  branch.  The  signal  beam  and  the 
reference  beam  combine  at  the  bearrKombiner  and  are  Fourier  transformed  onto  the  spatial  light 
modulators,  which  filter  the  combined  beams,  and  then  arc  collected  onto  an  area  photodetector 
The  area  photodetector  collects  all  light  from  the  signal  beam  and  reference  beam.  The 
photodetector  integrates  tire  li^t  intensity  over  its  surface  producing  two  bias  currents,  gif(t)  and 
g2f(t),  and  the  desired  cross-product  current  g3f(t). 

Figure  3.2  shows  the  same  optical  signal  processor  but  with  the  light  detected  in  tire  image 
plane  of  the  acousto-optic  cell  A  point  photodetector  that  is  just  large  enough  to  accept  the 
reference  beam  is  used  to  detea  the  light;  the  photodetector  does  not,  however,  accept  all  of  the 
signal  beam  because  it  is  spread  over  the  entire  image  of  the  acousto-optic  cell.  The  intensity  of  the 
light  is  integrated  over  the  surface  of  the  photodetector  to  give  the  two  bias  currents,  g^  i(t)  and 
g2i(t),  and  the  aoss-product  current  g3i(t).  It  rqjpears  that  g3{<t)  is  greater  than  g3i(t),  since  less 


20 


lig^t  is  collected  by  the  point  photodetector  in  die  image  plane.  We  give  an  exanqile  of  a  cw 

interferer  to  illustrate  die  apparent  discrepancy. 

From  Section  32  we  recall  that,  whoi  the  acousto-optic  cell  is  uniformly  illuminated,  the 

Fourier  transform  of  a  cw  signal  is  a  sine  distribution.  To  successfully  notch  die  cw  interferer  we 

remove  at  least  the  mainlobe  of  the  sine  function.  Hgure  3.6  shows  the  sine  function  after  the 

mainlobe  is  successfuHy  notched,  with  the  area  photodetector  collecting  the  energy  that  is  not 

blocked  by  the  spatial  light  modulators.  We  now  Fourier  transform  die  sidelobe  energy  and  detect 

the  signal  in  the  image  plane  of  the  acousto-optic  cell  using  a  point  photodetecton  Figure  3.7 

shows  the  inverse  Fourier  transform  of  the  sidelobes  of  the  sine  function  with  the  corresponding 

point  photodetector  collecting  the  light  The  point  photodetector  can  therefore  be  about  1/200  the 

size  of  the  signal  beam  in  the  image  plane,  since  it  only  needs  to  be  as  large  as  the  resolution  size  of 

the  system,  and  the  acousto-optic  cell  used  in  the  Presort  Processor  has  a  time  bandwidth  product 

of  200.  We  observe  that  only  a  small  portion  of  the  signal  beam  in  the  image  plane  is  detected  by 

die  point  photodetector.  Therefore,  it  appears  that  the  photocurrents,  and  hence  the  notch  depths, 

are  different  when  detecting  the  light  in  the  Fourier  plane  as  opposed  to  the  image  plane.  We  now 

show  that  the  cross-produa  photocurrents  are  equal  and  explain  the  apparent  discrepancy. 

3J.1  Derivation  of  the  photocurrent  using  an  area  detector  in  the  Fourier  plane 

We  derive  the  photocurrents  when  Fourier  plane  detection  is  used.  We  represent  the  signal 

directly  after  the  acousto-optic  cell  as  f(x,t),  and  the  reference  beam  as  r(x).  We  Fourier  transform 

both  beams,  and  represent  their  Fourier  transforms  as  F(a,t)  and  R(a).  We  represent  the 

amplitude  transmittance  of  the  spatial  light  modulator  as  H(a).  After  being  modified  by  the  spatial 

light  modulator,  the  signal  beam  is  given  by 

F(aj)H(a). 

Similarly,  the  modified  reference  beam  is  represented  as 

R(a)H(a). 

The  intensity  of  the  light  incident  on  tiie  photodetector  is 

/(a,r)  =  |F(a,/)//(a) + R(a)//(af. 


(3.28) 
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We  expand  equation  (3^8)  to  get 

l(a.l)  =  |F(o,()«{of  +  ^la)H(af + 2  Re{f  (o,/)i/(a)«*(a)«'(o)}.  ^9) 

We  integrate  die  intensity  over  the  photodetector  surface  to  give 

(3  30) 

where  we  assume  a  uniform  response  over  the  surface  of  the  photodetector,  P(a)  equals  one  in  the 
range  of  a(.-aco^a<iac+aco,  and  is  equal  to  zero  elsewhere.  We  also  assume  that  the  signal  and 
the  spatial  light  modulator  amplitude  transmit^ce  are  bandlimited  to  ac-aco^0i^c'''^co' 


We  integrate  die  signal  bias  term  over  the  surface  of  the  photodetector  to  get 

ft/ W  =  r‘*‘“nc‘.ma)F’(a,t)H\a)da. 


(3.31) 


or 


(3.32) 


where  we  relax  the  limits  of  integration  to  infinity  since  the  signals  of  interest  are  bandlimited  to 
ac-aco^a:Sag+a5o.  In  the  following  derivations  throughout  the  chapter,  we  assume  that  the 
signals  are  bandlimited  so  the  photodetectorcoUects  ail  the  light,  unless  otherwise  stated. 
Similarly,  we  integrate  the  reference  bias  term  over  the  surface  of  the  photodetector  to  give 

(3.33) 

where  the  limits  on  the  integration  are  infinite  since  all  the  reference  beam  and  signal  beam  light  is 
collected.  We  integrate  the  cross-product  term  to  get 

«3/(»)  =  J^2Re{f(o,;)H(a)«‘(o)H‘(a)]da. 

or 

= £.2Re{f  (<*.')«’(a)|ff(of  }<«•  (3  35) 

We  have  corrqileted  the  derivation  of  the  photocunents  when  Fourier  plane  detection  is  used.  We 

now  show  that  the  photocunents  represented  by  equations  (3.32),  (3.33),  and  (3.3S)  are  equal  to 

the  photocunents  when  area  detection  is  used  in  the  image  plane. 

33^  Derivation  of  the  photocurrents  using  an  area  detector  in  the  image  plane 
We  perform  a  Fourier  transform  on  the  filtered  beams  and  detea  the  light  in  the  image 
plane.  The  reference  beam  in  die  image  plane  is  represented  by 

(3.36) 


J2^R(a)H(ay^’^da, 
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where  we  perform  the  Fourier  transform  with  respect  to  die  firequency  variable  a.  The  signal  beam 
is  represented  by 

We  integrate  the  intensity  over  the  surface  of  the  area  photodetector  to  get  the  signal  bias  term  as 

(3.38) 

where  we  have  introduced  a  dummy  variable  y  in  the  second  integral.  We  rearrange  die  integrals  to 


=  (3  39) 

We  integrate  over  u  to  get 

=  2Re|j^j^F(a,r)//(a)F*(y.r)//*(r)5(a-  y)rf05dy}. 

We  use  the  sifting  property  of  the  delta  function  to  get 

«ii(')=|l|f(r.f)«(r)|V  (3  4J) 

We  observe  that  the  signal  bias  term,  when  image  plane  detection  is  used,  is  equal  to  die  signal  bias 
term  when  Fourier  plane  detection  is  used,  since  equation  (3.41)  equals  equation  (3.32). 


(3.41) 


We  use  the  same  integration  technique  to  find  the  reference  beam  bias  term: 

82M^Zy^{r)tf{ridy- 


(3.42) 


Equation  (3.42)  is  equal  to  equation  (3.33),  therefore  showing  that,  when  area  detection  is  used  in 
the  Fourier  plane  and  area  detection  is  used  in  the  image  plane  to  detect  the  light,  the  bias  terms  are 


equal. 


We  now  consider  the  cross-product  term  which  is  given  by 
g3i(0  =  j^^2Re|j7^F(a,r)//(a)e^^*“da|^/?*(y)//*(y)e"-'^*’*‘dy|j£i«. 


We  rearrange  die  order  of  integration  in  equation  (3.43)  to  get 
«3,(t)  =  2Re|  f  r„  r.  F(g.')H(g)«‘(r)tf 
We  perform  the  integration  over  u  to  get 

g3,(r)  =  2Re|£‘^  J^F(a,r)//(a)F*(y)//*(y)(5(a  -  y)dady|. 


(3.43) 


(3.44) 


(3.45) 
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We  now  use  the  sifting  property  of  the  delta  function  as  we  integrate  over  a  to  get 

Equation  (3.46)^uals  equations  (3.35),  we  therefore  have  shown  that  if  all  the  light  is  collected, 

frequency  plane  detection  and  image  plane  detection  are  equivalent 

3.4  Derivation  of  the  photocurrents  using  a  point  detector  in  the  image  plane 

In  the  previous  derivations  we  relaxed  tte  limits  of  integration  to  infinity  since  all  the  light 
is  collected.  We  recall  that  a  point  photodetector  in  the  image  plane  is  sufficient  to  collect  all  the 
light  for  which  tire  reference  and  signal  beams  overlap.  The  point  photodetector  in  the  inu^e  plane 
does  not,  however,  accqrt  all  of  tire  signal  beam,  and  we  carmot  set  the  limits  of  integration  to 
infinity.  In  the  following  derivation,  we  begin  by  letting  the  photodetector  have  a  width  of  L, 
where  L  is  large  enough  to  collect  all  the  reference  beam  and  signal  beam;  L  is  the  length  of  the 
acousto-optic  cell.  The  cross-product  photocunent  in  tiie  image  plane  is  given  by 
g2i{t)  -  J^t^2Re|jJ^f 

(3.47) 

where  we  represent  the  first  integral  by  f  (u,i)  and  tire  second  Fourier  integral  by  r'(u).  The  cross- 
product  photocurrent  is  then  represented  as 


g3,(f)  =  2Re|j^/(u,r)r'(u)du|, 


(3.48) 


f  {u,ty  iu)duj 

where  equation  (3.48)  is  equal  to  equation  (3.35),  since  L  accepts  all  the  light  in  the  reference  and 
signal  beams.  We  recall  that,  by  defirution,  we  obtain  a  cross-product  photocurrent  where  tire 
signal  beam  and  reference  beam  overlap.  Since  the  reference  beam  is  much  smaller  than  the  signal 
beam,  we  can  let  the  size  of  the  photodetector  shrink  to  W,  where  W  equals  the  extent  of  the 
reference  beam  and  W«L.  The  cross-product  photocunent  then  becomes 

g3, (f)  =  2  Rej f  iu^ty  (u)du|. 

We  therefore  discover  tirat  the  cross-produa  photocunent  represented  by  equation  (3.49)  is  equal 
to  equation  (3.48),  which  is  equal,  in  turn  to  equation  (3.35),  where  equation  (3.35)  represents  the 
photocurrent  when  aU  the  light  is  collected  in  the  Fourier  plane  using  an  area  photodetector  We 


(3.49) 
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.  have  therefore  shown  Aat  the  cross-product  photocurrents  are  equal  whether  we  use  a  point 
-  photodetector  in  the  ima^  plane  or  an  area  photodetector  in  the  Fourier  plane.  There  is  no 
fundamental  difiercnce,  as  far  as  die  cross-product  tenn  is  concerned,  between  Fourier  plane 
detection  and  image  plane  detection,  provided  that  all  the  overh^ing  and  colinear  light  is  collected. 
However,  since  we  clearly  do  not  collect  all  the  signal  beam  energy  when  point  detection  is  used. 


what  affect  does  diis  have  on  the  system  performance? 


The  signal  bias  tenn,  when  detected  with  a  finite  width  photodetector  L  is  given  by 


(3.50) 


We  represent  the  first  Fourier  integral  as  f  (u,t)  and  the  second  as  its  complex  conjugate. 
Therefore,  the  signal  bias  term  photocurrent  becomes 


'  (3.51) 

We  now  let  the  photodetector  size  shrink  to  W  so  that  it  is  just  large  enough  to  collect  all  the 

reference  beam.  Equation  (3.51)  is  therefore  represented  as 


Since  the  photodetector  is  smaller  than  the  extent  of  the  signal  beam  we  get 

.w  .  2 


I  small  portion  of  the  signal  beam, 

W  1  fit  ,  ^ 

J\|/'(u,/)|  du<p,[/’(M,r)|  du 


Since  we  are  integrating  over  a  small  portion  of  the  signal  beam,  we  find  that 

w  1  L  2 


(3.52) 


(3.53) 


(3.54) 


where  we  assumed  that  W«L.  We  therefore  observe  that  dw  signal  beam  bias  tom  is  reduced 
when  a  point  photodetector  is  used  in  the  image  plane  to  detect  die  light  As  a  consequence, 
collecting  a  small  portion  of  the  signal  beam  reduces  die  shot-noise  of  the  optical  system,  therefore 
improving  the  signal-to-noise  ratio,  as  we  show  in  the  next  section. 

Following  a  similar  argument  as  presented  above,  and  recalling  that  the  point  photodetector 
is  large  enough  to  collea  all  the  reference  beam,  we  find  that  the  reference  beam  bias  term  is 

.L  .  .2  w  .  .2 


^2iW  =  |t|r’(M,t)(  du=J\|r'(u,f)|  du. 


(3.55) 
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Since  all  the  reference  beam  is  collected  in  both  the  Fourier  plane  using  an  area  photodetector  and 
the  image  plane  using  a  point  photodetector,  we  find  that  the  reference  beam  bias  terms  are  equal 
under  the  two  measurement  conditions. 

3.5  The  effect  of  the  detection  schemes  on  the  signal>to-noise  ratio 


The  signal-to-noise  ratio  at  the  photodetector  ouq)ut  is  given  by 


SNR^ 


2eB(i^  +  i^)Ri  +  4kTB  ’  (3  5g) 

where  i,  is  die  signal  current,  e  is  the  charge  on  an  electron,  B  is  the  post-detection  bandwidth,  i^  is 


the  dark  current,  i,^  is  the  average  signal  current,  k  is  Boltzmann’s  constant,  T  is  die  temperature  in 
degrees  Kelvin,  and  is  die  photodetector  resistance.  The  average  signal  current,  i,.,  is 
proportional  to  die  power  of  the  reference  beam  plus  the  power  of  the  signal  beam  plus  the  power 
of  the  cross-produa  term.  Since  the  cross-product  term  is  constant  for  these  calculations,  we 
examine  just  the  contribution  from  the  bias  terms. 

We  present  an  example  of  the  effect  on  the  signal-noise  ratio  to  facilitaie  our  discussion  of 
the  two  detection  schemes.  Figures  3.1  and  3.2  show  an  optical  processor  where  the  light  is 
detected  in  the  Fourier  and  image  planes  of  the  acousto-optic  cell,  respectively.  We  assume  the 
beamsplitters  are  50/50  power  splitters  and  that  the  laser  has  an  ouqiut  power  of  10  mW.  In  the 
Fourier  plane  all  the  reference  and  signal  beams  are  collected  by  an  area  detector  For  200 
simultaneous  cw  inputs  to  die  acousto-optic  cell  widi  a  di&action  ratio  of  one  microwatt  of  optical 
power  per  frequency,  we  get  200  microwatts  of  signal  power  at  the  Fourier  plane;  these  optical 
power  values  arc  from  measurements  made  on  the  Presort  Processor  We  have  a  total  referrace 
beam  power  of  2.5  mW  in  the  Fourier  plane.  The  total  power  incident  on  the  area  photodetector  in 
the  Fourier  plane  is  therefore  2.5  mW  +  0.2  mW  =  2.70  mW. 

At  the  image  plane  we  collect  the  reference  beam  plus  a  portion  of  die  signal  beam;  we 
assume  that  the  photodetector  is  1/200  the  size  of  die  signal  beam.  The  total  power  incident  on  the 
point  detector  is  then 


2.50  mW  200uW/200  =  2.501  mW. 
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We  now  calculate  the  percentage  power,  relative  to  the  area  photodetector  power,  diat  is  incident 
'Upon  the  point  photodetecton  The  percentage  power  is 


2.501 

2.700 


X 100  *92.6%, 


(3.57) 


in^lying  diat  point  detection  does  not  have  a  significant  advantage  over  Fourier  plane  detection  in 
terms  of  shot  noise  reduction.  We  note  diat  die  type  of  detecticHi  scheme  used  is  dqiendent  on  the 
type  of  signals  that  are  to  be  processed.  We  further  note  that  using  point  detection  over  area 
detection  would  be  advantageous  when  die  system  is  shot-noise  limited  and  the  signal  and 
reference  beams  are  the  same  strength  at  die  photodetectOL  Finally,  point  detection  is  better 
because  its  c^acitance  is  smaller  than  the  capacitance  for  area  detectors,  and  it  can  tlwrefore  be 
used  for  wider  bandwidth  signals. 

3.6  Discussion  and  conclusions 

Figure  3.7  shows  an  envelope  of  a  cw  signal,  at  the  image  plane,  when  the  mainlobe  of  the 
cw  interferer  is  excised  in  the  frequency  plane.  The  small  sine  function  in  the  middle  of  the  graph 
is  the  image  plane  representation  of  the  reference  beam;  the  reference  beam  and  the  signal  beam 
have  an  amplitude  of  one.  Recall  that  the  cross-prcxiuct  photocunent,  g3i(t),  is  proportional  to  the 
integral  of  the  reference  beam  multiplied  by  the  signal  beam.  The  reference  beam  is  positioned  on 
the  optical  axis  (i.e.,  xo=0).  If  we  were  to  slide  the  reference  beam  to  the  left,  across  the  image 
plane,  we  observe  diat  the  reference  beam  multiplied  by  the  signal  beam  decreases  until  a  minimum 
value  is  reached  at  the  zero  point;  the  minimum  value  irrplies  that  the  notch  depth  is  a  maximum.  If 
we  normalize  the  image  plane  to  a  length  of  one  with  the  outer  limits  being  0.50  and  -0.50,  the 
zero-crossing  occurs  at  plus  or  minus  0.215  from  the  optical  axis.  This  result  agrees  with 
Erickson’s  work  [1],  in  which  he  showed  that,  for  the  special  case  of  rectangular  apodization  with 
three-fourths  of  die  mainlobe  excised,  a  deep,  narrow  notch  occurs  at  0. 17.  Since  the  entire 
mainlobe  is  excised  in  our  example,  the  distance  of  the  zero-crossings  from  die  optical  axis  are 
increased.  By  using  this  image  plane  representation  we  therefore  gain  some  insight  into  the 
relationship  between  the  notch  depth  and  the  reference  beam  position. 
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Figure  3.8  shows  the  image  plane  when  the  mainlobe  and  neighboring  sidelobes  are 
excised;  plot  #1  is  when  the  mainlobe  is  excised,  plot  #2  is  when  the  mainlobe  and  the  nearest 
sidelobes  are  excised,  and  plot  #3  shows  when  the  mainlobe  and  the  nearest  two  sidelobes 
excised.  We  observe  diat  die  photocmrent  decreases  as  successive  sidelobes  are  excised  if  the 
reference  beam  is  on  the  optical  axis.  We  hirther  note  that  there  are  a  variety  of  zero-points  for  the 
reference  beam  to  dwell  as  more  sidelobes  are  excised. 

We  now  explore  the  performance  of  an  optical  excisor,  namely  the  Presort  Processor, 
which  uses  frequency  plane  notch  filtering  to  excise  narrowband  interference  fit>m  wideband 
signals.  We  then  design  and  test  a  system  diat  modifies  die  Presort  Processor  so  that  the  frequency 
excision  is  open-loop  adaptive. 
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4.  Presort  Processor  modification 

In  this  chapter  we  discuss  the  operadon  of  the  Presort  Processor  in  detail;  this  includes  the 
characteristics  of  the  signal  acousto-optic  cell,  the  spatial  light  modulators,  and  a  description  of  the 
frequency  excision  technique  used  We  describe  how  we  modified  die  Presort  Processor  for  diis 
study;  an  analysis  is  given  of  die  effect,  on  the  Presort  Processor,  of  inserting  an  optical  window 
in  the  interferometer  to  split  the  signal  beam.  We  describe  the  operation  of  die  post-processing 
electronics  which  are  used  to  determine  the  location  of  the  narrowband  interferers. 

4.1  D<«cription  of  the  Presort  Processor 

The  Presort  Processor  is  an  optical  signal  processor,  developed  by  the  Harris  Corporation, 
to  excise  narrowband  interference  from  wideband  signals;  the  description  of  the  Preson  Processor 
in  Sections  4.1  through  4.12  is  based  upon  information  obtained  from  the  Final  Report  [15]. 
Figure  4.1  shows  the  basic  components  of  the  Presort  Processor.  The  solid  state  laser  has  a 
wavelength  of  830  nm,  with  a  maximum  power  ouqiut  of  40  mW.  The  laser  includes  collimating 
optics  which,  toge  Jier  with  lenses  and  L2,  expand  the  beam  to  a  4  mm  length  in  the  direction 
parallel  to  the  page.  The  beam  is  250  |im  wide  in  the  direction  perpe  dicular  to  the  page. 

The  Presort  Processor  is  based  on  a  Mach-Zehnder  interferometer  with  an  acousto-optic 
cell  and  a  prism  in  the  lower  branch  of  the  interferometer,  and  with  a  lens  in  the  upper  branch.  The 
prism  expands  the  light  in  the  direction  of  acoustic  propagation  from  a  4  mm  beam  to  an  8  mm 
beam  to  illuminate  the  acousto-optic  cclL  The  acousto-optic  cell  receives  the  electrical  signal,  f(t), 
which  is  distorted  by  narrowband  interferers,  and  modulates  tire  light  in  both  space  and  time. 

The  light  in  the  upper  branch  of  the  interferometer  passes  throu^  lens  L3  to  create  a  point 
source,  12  pm  in  size,  i^ch  is  located  exactly  the  same  distance  from  tiie  beam  combiner  as  is  the 
acousto-optic  cell.  The  point  source  acts  as  a  probe  on  the  acousto-optic  cell,  causing  the 
wideband  photodetector  to  detea  that  part  of  the  signal  in  the  acousto-optic  cell  which  overlaps  the 
point  source.  To  further  explain  this  operation,  we  approximate  the  point  source  as  a  delta 
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function.  At  the  photodetector,  which  is  in  an  image  plane  of  the  acousto-optic  cell,  we  have  the 
convolution  of  the  signal  with  a  delta  function,  which  recovers  the  signal,  as  we  discussed  in 
Section  3.2.1.  — 

Since  both  segments  after  die  second  beamsplitter  operate  identically,  we  describe  just  the 
strai^t  segment  in  detail.  Spherical  lenses  L4  and  display  die  Fourier  transform  of  the  signal 
f(t)  at  the  plane  of  the  spatial  light  modulator.  The  spatial  light  modulator  is  used  to  block  or  pass 
various  parts  of  the  spectrum  of  the  signal.  The  spadal  light  modulators  used  in  this  system  have  a 
50  percent  duty  cycle,  resulting  in  half  the  spectrum  being  covered  by  SLMl  in  the  straight 
segment,  while  the  odier  half  of  the  spectrum  is  covoed  by  SLM2  in  the  folded  segment 
Cylindrical  lens  L5  increases  the  waist  of  the  beam  from  250  |im,  at  the  acousto-optic  cell,  to  400 
|im  at  the  spatial  light  modulator  After  the  beam  passes  throu^  the  spatial  light  modulator  it  is 
collected  by  spherical  lens  L7,  which  focuses  the  light  onto  a  wideband  photodetector.  The 
photodeteaor  output  from  the  straight  segment  is  then  added  to  the  photodetector  output  of  the 
folded  segment,  so  that,  in  effect,  the  two  spatial  light  modulators  cover  the  entire  bandwidth. 

4.1.1  Te02  longitudinal  phased  array  acousto-optic  cell 

The  signal  acousto-optic  cell  in  the  lower  branch  of  the  interferometer  is  a  Te02  longitudinal 
mode,  phased  array  device.  Figure  4.2  shows  the  passband  characteristics  of  the  200  MHz 
bandwidth  cell;  the  bandwidth  of  the  cell  ranges  from  250  MHz  to  450  MHz.  The  time-bandwidth 
product  of  the  cell,  given  that  the  length  of  the  cell  is  8  mm,  that  the  acoustic  velocity  of  the  cell  is 
4200  m/s,  and  tiiat  the  bandwidth  is  200  MHz,  is  340.  The  cell  is  illuminated  by  a  Gaussian 
distribution  which  reduces  the  number  of  resolvable  friequencies  by  a  factor  of  approximately  1.5 
because  the  Fourier  transform  of  a  truncated  Gaussian  illumination  has  a  wider  mainlobe  than  that 
for  uniform  illumination.  The  number  of  resolvable  frequencies  is  therefore  200,  which  is  equal  to 
the  numhrr  of  spatial  li^t  modulator  elements. 
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4.1^  Spatial  light  modulators 

Figure  4.3  shows  the  spatial  light  modulators,  which  consist  of  two  100-chaonel  Te02 
longitudinal  mode,  acousto-optic  cells.  The  transducers  for  each  ceU  are  on  center-to-center 
spacings  of  250  |xm,  with  transducer  widths  of  125  pm,  and  a  ceU  height  of  10  mm.  These  spatial 
light  modulators  have  a  50  percent  duty  cycle  resulting  in  both  segments  of  the  Presort  Processor 
'  being  dedicated  to  the  fipequency  excision.  The  spatial  light  modulator  resolution  in  the  fiiequency 
plane  is  equal  to  the  width  of  one  transducer,  which  is  do=1254m.  Figure  4.3  shows  how  the 
multi-channel  acousto-optic  cells  are  interleaved  so  that  one  ^atial  light  modulator  covers  the  odd 
frequency  bins  and  the  other  covers  the  even  frequency  bins.  If  we  were  to  slide  me  spatial  light 
modulators  vertically  along  the  dashed  lines  undi  they  overlap,  every  frequency  bin  between 
250  MHz  and  449  MHz,  with  a  resolution  of  1  MHz,  is  represented  by  a  spatial  light  modulator 
element  If  a  narrowband  interferer  is  present  the  appropriate  acousto-opdc  channels  are  turned 
ofr,  and  the  light  which  corresponds  to  that  interferer  is  no  longer  incident  on  the  photodetectors. 

Hgure  4.4  shows  the  theoretical  notch  shape  provided  by  the  spatial  light  modulators  to 
notch  a  narrowband  interferer,  directly  below  the  grs^h  are  the  locations  of  the  spatial  light 
modulator  elements.  Curves  one  and  two  show  the  notch  shape  with  one  and  two  spatial  light 
modulator  elements  swimhed,  respectively.  Notice  that  switching  three  consecutive  spatial  light 
modulator  elements  gives  an  excision  dq>th  of  25  dB,  while  switching  five  elements  provides  a 
.notch  depth  of  28  to  33  dB. 

We  use  an  example  to  illustrate  some  of  the  characteristics  of  the  spatial  li^t  modulators. 
We  assume  mat  a  25  dB,  referenced  to  me  noise  level,  1  MHz  bandwidth,  narrowband  interferer  is 
distorting  a  spectrum.  To  excise  me  25  dB  narrowband  interferer  we  need  to  switch  at  least  three 
spatial  light  modulator  elements  as  Figure  4.4  shows;  if  a  N  MHz  wide  notch  is  desired,  N+1 
spatial  light  modulator  elements  must  be  switched.  Figure  4.4  shows  that  switching  two  spatial 
light  modulator  elements  results  in  a  22  dB  notch  and  is  merefoie  not  sufficient  to  perform  me 
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excision.  Therefore  to  excise  a  25  dB  interferer,  at  least  three  spatial  light  modulator  elements  must 

be  switched. 

4.2  Modification  of  the  Presort  Processor 

The  Presort  Processor  currently  has  a  panel  of  200  switches  to  control  the  spatial  light 
modulators.  We  are  implementing  a  power  spectrum  analyzer,  along  with  post-detection 
electronics,  to  provide  the  control  signals  needed  to  adaptively  switch  the  spatial  light  modulator 
elements.  Figure  4.5  shows  how  the  Presort  Processor  is  modified  to  split  the  signal  beam  to 
implement  a  power  spectrum  analyzer  The  an:q)litude  of  the  signal  beam  is  divided  by  an  optical 
window,  which  is  inserted  at  a  twenty  degree  angle  after  the  acousto-optic  cell,  to  reflect  4  percent, 
while  transmitting  96  percent,  of  the  light  The  reflected  signal  beam  is  then  Fourier  transfonned 
in  the  direction  of  propagation  of  the  acousto-opdc  cell,  and  is  imaged  in  the  perpendicular 
direction  by  lenses  Lj  2.  Lj  3,  and  L14.  The  spectrum  is  incident  on  a  CCD  photodetector  array,  and 
sequentially  clocked  into  the  post-detection  electronics.  The  description  of  the  operation  of  ibe 
post-detection  electronics  is  deferred  until  Section  4.5.1. 

When  we  split  the  signal  beam  with  the  optical  window,  the  beam  passing  through  the 
optical  window  is  displaced  in  the  direction  of  propagation  of  the  acousto-optic  celL  Furthermore, 
the  second  surface  reflection  from  the  optical  window  produces  a  ^ost  beam  parallel  to  the  first 
surface  reflection,  which  causes  a  fringe  pattern  to  occur  on  the  spectrum.  Finally,  the  wedge  on 
the  optical  window,  caused  by  die  nonparallel  surfaces  of  the  optical  window,  causes  die  spectrum 
of  the  signal  to  shift  at  the  spatial  light  modulators.  These  effects  are  explored  in  the  following 
sections. 

4,2.1  Calculation  of  the  required  clear  aperture 

In  this  section  we  calculate  the  window  clear  aperture  needed,  based  on  the  length  of  the 
optical  beam  after  the  acousto-opdc  cell  The  illuminated  part  of  the  acousto-optic  cell  has  a  length 
of  8  mm,  and  the  beam  is  diffracted  by  the  acousto-optic  cell  at  a  maximum  angle  of  6|M).0198 
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radians  with  respect  to  the  center  frequency  (350  MHz).  Figure  4.6  shows  how  the  optical 
window  is  inserted  at  a  20  degree  angle  with  respect  to  the  optical  axis.  The  required  height  of  the 
beam  at  the  optical  window,  which  is  approximately  100  mm  from  the  acousto-optic  cell,  is  the 
sum  of  tile  apertures  needed  to  collect  the  diffracted  and  undiffracted  light: 

ho=2(100  mm)(0.0198  radians)-i-8  miir=11.95  mm. 


Hence,  the  clear  aperture  needed  for  this  configuration  is 

lL95mm 

CA  - - -  =  12.72inm. 

sin(90  -  20)°  (4.1) 

The  clear  aperture  for  an  optical  window  of  length  30  mm,  when  inserted  at  a  20  degree  angle,  is 

21.75  mm.  Since  the  clear  tqierture  of  the  optical  window  is  larger  tiian  the  hei^t  of  the  optical 

beam,  a  30  mm  diameter  optical  window  is  sufficient 

4.2.2  Displacement  of  the  second  surface  reflected  wave 

Figure  4.6  shows  the  angle  of  incidence  6|,  which  is  twenty  degrees,  on  the  optical 

window.  We  know  from  Snell’s  law  that 

«isin0i  =sn2sin02.  (4.2) 

where  n^  and  n2  are  the  indices  of  refraction  of  the  first  medium  and  second  medium,  respectively, 

and  0 1  and  62  are  the  angle  of  incidence  and  the  angle  the  transmitted  ray  makes  with  respect  to  the 

normal,  respectively.  In  our  case  n^  is  equal  to  one,  and  02  is  equal  to  the  index  of  refraction  of  the 

fused  silica,  which  is  1.45. 


We  first  solve  for  62  in  equation  (4.2)  to  get 


n. 


sin  ft 


^"2 


(4.3) 


$2  5=  sin~*‘ 

Since  the  thickness  of  the  optical  window  is  equal  to  t,  the  conq>onent  parallel  to  the  optical 
window  that  the  beam  travels  is 

xi=rtana2.  (4.4) 

At  the  second  surface,  part  of  the  light  is  transmitted,  pan  of  the  light  is  reflected,  and  the 


remaining  light  is  absorbed  by  the  optical  coating.  From  the  law  of  reflection,  we  know  that  the 
angle  of  incidence  is  equal  to  the  angle  of  reflectance,  therefore,  the  total  distance  the  beam  travels 
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along  die  optical  window  is  equal  to  2x2. 

The  reflected  beam  from  die  back  surface  now  escapes  from  the  optical  window  and  makes 
a  20  degree  an^  widi  respect  to  die  normal  of  the  front  surface.  Hence,  the  ray  reflected  from  the 
first  surface  is  parallel  to  the  ray  which  escapes  from  the  optical  window  from  the  second  surface 
reflection;  for  now  we  are  assuming  the  two  surfaces  of  die  optical  window  are  parallel 
Therefore,  die  angle  the  beam  makes  with  die  optical  window  is  70  degrees,  and  the  total 
displacement  of  the  beam,  X3,  can  be  calculated  by  noting  that 

X3=X2COSdi.  (4  5) 


JC3  s2/tan| 


• 

sin“^ 

-^sin0i 

. 

J. 

COS  01, 


(4.6) 


where  the  total  displacement  is  in  terms  of  the  input  angle,  die  thickness  of  die  optical  window,  and 
the  indices  of  refraction.  Since  the  optical  window  is  6  mm  thick,  we  substitute  the  parameters  into 
equation  (4.6)  and  find  that  the  total  displacement  is  X3»2.74  mm. 

We  now  find  the  total  displacement  of  the  transmitted  beam,  and  we  denote  this  quantity  as 
X5.  We  see  that  the  optical  path  the  beam  traverses  throu^  the  material  is 


t 

^4  - - TT- 

COS  92 

Furthermore,  we  recognize  that  the  displacement  of  the  transmitted  beam  is 

a:5=X4sin(0i-02)- 


(4.7) 

(4.8) 


We  substitute  equations  (4.2)  and  (4.7)  into  equation  (4.8)  and  we  get 

#  r  ^  ^ 

t  .  Irt  _.  -1 


X5=- 


•  -1  •  n 

cosi  sm  — *-sm0i 


smi 


01  -  sin 


\n2 


'  sin  01 


(4.9) 


We  substitute  the  values  of  the  paraoKters  into  equation  (4.9)  to  find  that  X5=0.684  miiL 

In  summary,  we  find  that  the  second  surface  reflected  wave  is  displaced  by  2.74  mm  with 
respect  to  the  first  surface  reflected  wave,  and  that  the  transmitted  beam  is  displaced  by  0.684  mm 
with  respea  to  the  incident  beam  on  die  optical  window.  We  now  explore  the  effect  these 
displacements  have  on  the  powa  spectrum  at  the  CCD  photodetector  array. 
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4.2.3  Analysis  of  the  second  surface  reflection 


We  now  model  the  first  and  second  surface  reflections  of  die  optical  window.  We  let  the 
beams  of  li^t  reflected  by  the  optical  window  be  represented  by  g(x).  Therefore,  the  reflected 
light,  g(x),  is  represented  by 

gix)  =  gi(x)+g2M,  (4  10) 

where  g^fx)  and  g2(x)  are  the  first  and  second  surface  reflections,  respectively.  The  second 

surface  reflection  has  the  same  content  as  the  first  surface  reflection  but  is  physically  displaced,  has 

an  angular  shift,  and  is  attenuated  with  respect  to  the  first  surface  reflection.  The  second  surface 


reflection  can  therefore  be  represented  as  a  ftmcdon  of  the  first  surface  reflection: 

82 (X)  =  bgi(x - (4  j 

where  b  is  the  attenuation  constant,  X3  is  the  physical  displacement,  and  equals  Q^/X.  We 

substitute  equation  (4.1 1)  into  equation  (4.10)  to  get 

8(x)  =  81  (x) + bgi  (X  -  X2  .  (4  1 

The  function  g(x)  now  passes  through  a  scries  of  lenses  which  perform  the  Fourier  transform. 


The  Fourier  transform  of  g(x)  is 

0(a)  =  G,(oi)+iiG,(a-a,y^”“>.  (4 

where  the  exponential  term  is  caused  by  the  physical  displacement  of  the  two  reflected  beams. 


The  incident  light  falls  on  a  photodcteaor  array,  which  detects  the  intensity  of  G(a): 
m  =  |G,(a)f  +  i.2|Gi(a  -  +2Re(Gi(a)K;i‘(a  -  ).  ^4  ,4^ 

The  desired  spectrum  is  Gi(a)Gi*(a),  while  die  other  tvo  terms  are  unwanted. 

First  we  analyze  the  effect  of  the  physical  shift,  ^.a^F,  of  the  ghost  Fourier  transform 
which  is  caused  by  the  second  surface  reflection.  We  calculate  the  angular  deviation  from  the 
specifications  given  in  the  Melles  Griot  catalog.  The  02  WLQ  306  optical  window  has  a 
paraUelism  of  five  arcseconds;  five  arcseconds  equals  5/3600  degrees  or  2.42(10)"^  radians.  Since 
the  second  surface  reflection  occurs  on  a  tilted  surface,  the  angular  shift  with  respect  to  the  first 
surface  reflection  is  twice  the  angular  error  of  the  two  surfaces,  so  that  the  total  angular  shift  in 
radians  is 


06=2(2.42)(1O)-5  *  4.85(10)*5radians. 
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We  now  determine  the  effect  of  65  in  die  Fourier  plane  using  the  transfonnation 

4^  =  OfiAF  =  AF  =  dfiF  =  4.8482(10r^263.1mm)  =  \2.16taxi, 
where  the  focal  length  F=263. 1  is  calculated  in  Section  4.5.  The  angular  tilt  gives  a  physical  shift 
of  12.76  p.m  which  is  qiproximately  the  size  of  a  CXD  photodetector  element  (13  pm)  or  one- 
fourth  of  a  fircquency  resolution  element 

The  effect  of  die  physical  displacement  of  die  two  reflected  beams  is  now  analyzed.  Since 


the  shift  is  tqiproximately  one-fourth  of  the  frequency  resolution  size,  we  let 

Gi(a  -  a^) »  Glia)  (4.15) 

which  simplifies  the  analysis.  We  use  diis  approximation  in  equation  (4.14)  to  get 

m = |G,(of + b%(af  +2Re(C,(o)iiG,'(ay^'“> ).  (4  jg) 

We  now  take  die  real  part  of  the  last  term  of  equation  (4. 16)  and  combine  like  terms  to  get 

/(a)  =  |Ci(a)f  [l + 6^  -»■  2hcos(2;rccc3)]. 

The  first  term  in  the  brackets  of  equation  (4. 17)  is  the  desired  terra  The  second  and  third  terms  are 


undesirable,  since  the  second  term  is  caused  by  the  second  surface  reflectioa  which  also  sets  up 


the  third  terra  The  magnitude  of  b  determines  if  tlx  third  term  is  of  major  concern.  If  b,  which  is 


the  normalized  second  surface  reflectioa  is  equal  to  one,  the  spectrum  of  the  signal  is  significandy 


distorted,  since  the  term  inside  the  brackets  is  zero  when  the  cosine  argumoit  equals  180  degrees. 


If  the  fiinge  pattern  is  not  resolvable  by  the  photodetector  array  the  third  term  in 
equation  (4.17;  doe .  not  affect  the  output  From  Section  4.22  we  calculated  X3  to  be  equal  to 


2.74  mm.  The  frequency  of  die  fringe  pattern  is  therefore 

2.74mm 


=  12.5Ab. 


^  AF  (830nmX263.1mm) 

Since  the  length  of  the  800  photodetector  elements  is  10.4  mm,  one  cycle  of  the  fiinge  pattern 

spans  6.2  photodetector  elements,  and  the  photodetector  array  easily  resolves  die  fringe  pattern. 

If  the  optical  window  is  placed  in  the  intaferometer  without  an  anti-reflection  coating  on 


the  second  surface,  the  reflection  is  large  compared  to  the  first  surface  reflectioa  We  coated  the 
second  surface  with  Hebbar  A)76  which  gives  a  reflectance  of  less  than  0.4  percent  for  Ught  whose 
incident  angle  ranges  fiom  zero  to  thirty  degrees.  The  830  nm  solid  state  laser  is  polarized  with  the 
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electric  field  perpendicular  to  the  acoustic  propagation  of  die  acousto-opdc  cell,  hence,  in  the 

classical  sense,  the  polarization  is  parallel  to  die  boundary  of  the  window.  To  calculate  die  amount 

of  light  leflected^y  the  optical  window  we  use  the  Fresnel  reflection  coefficient  [16] 

^  _  /i^cosfii  -n^cosg; 

^  /Ij  cos  ^  ^2  cos  @2  ^4  J 

where  ni  and  n2  are  the  indices  of  refraction  of  the  first  medium  and  second  medium,  respectively, 

and  6i  and  62  are  the  angle  of  incidence  and  die  angle  the  transmitted  ray  makes  widi  respea  to  die 

normal,  respectively.  In  our  case  6]  =  20  degrees,  62  =  13.64  degrees,  n^  -  1.0,  and 

02^1.45,  so  that 

tp  1=0.20. 

To  calculate  the  power  reflected,  we  square  tb^  reflectioo  coefficient  to  get 

Ri=  Tpi^  =0.040, 

which  is  approximately  4  percent  of  the  li^t  at  20  degrees  incidence.  The  power  reflected  by  the 
second  surface  is 

R2=rp22=(l-0.040)(0.004)(0.960)=0.0037. 

Therefore  we  have  diat  rp2  =0.061,  and  we  calculate  b  as 

h*^  =  ^^  =  0.305. 

Tpi  0.200 

We  return  to  equation  (4.17),  and  substitute  in  the  various  parameters,  to  find  that 

/(a)  =  1  +  0.0930  +  2(0.305)cos(2;rccc3) 


or 

/(a)  =  L0930  +  0.610cos(2;roct3).  (4  19) 

From  equation  (4.19)  we  find  that  the  maximum  value  of  the  fiinge  pattern,  relative  to  the 

mean,  is  1(a)  =  1.7  which  is  a  2.31  dB  increase.  The  minimum  value  of  the  intensity  pattern  is 

0.48  which  is  a  -3.16  dB  decrease  from  the  mean.  The  fringe  pattern  therefore  varies  by  5.5  dB 

which  shows  that  the  response  to  a  30  dB  narrowband  interferer  may  range  from  26.84  dB  to 

32.31  dB.  This  range  is  acceptable  because  we  are  considering  the  case  where  the  interferer  is 

20  dB  to  30  dB  above  the  signal  level,  and  we  do  not  need  fine  amplitude  resolution  to  discriminate 

between  the  signals. 

We  summarize  this  analysis  by  noting  the  presence  of  the  ghost  spectrum  at  the 
photodetector  array  distortion  to  the  original  spectruiiL  The  amount  of  distortion  is  dq)endent  on 
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the  strength  of  the  ^ost  spectrum,  and  on  the  physical  displacement  of  the  ^ost  spectrum  with 
respect  to  die  desired  qiectrum.  Since  the  ghost  spectrum  is  displaced  one-fourth  of  a 
resolution  elemmt  and  the  second  surface  reflecdon  is  attenuated  by  10.3  dB  with  respect  to  the 
first  surface  reflection,  the  amount  of  distortion  is  small.  The  fringe  pattern  on  the  spectrum  will 
cause  some  of  the  narrowband  interferers  to  appear  stronger  or  weaker,  however,  since  the  second 
surface  is  anti-reflection  coated,  this  effect  is  small. 

4.2.4  Effect  of  the  wedge  of  the  optical  window  at  the  spatial  light  modulators 
In  this  section  ive  explore  the  effect  of  the  wedge  on  the  optical  window  on  die  Fourier 

planes  at  the  spatial  light  modulators.  Recall  diat  die  optical  window  chosen  for  diis  experiment 
has  a  parallelism  of 

(t>  =  2.42(10)'^  rad. 

How  much  does  this  angular  tilt  affea  the  Fourier  plane?  We  know  that  9=Xa  so  that 

^  =  alF  =  0F.  (4.19) 

Hence,  the  amount  of  displacement  in  the  Fourier  plane  is 

A4  =  (634m/n)(2.42)(10)‘^  =  15.34Mm. 

Since  the  spatial  resolution  size  in  the  Fourier  plane  is  125  pm,  the  spectrum  is  shifted  by  about  1/8 
of  a  resolution  element  The  shift  can  be  easily  corrected  in  the  Presort  Processor  by  moving  the 
spatial  light  modulators  as  part  of  the  calibration  procedure. 

4.2.5  Analysis  of  a  shifted  signal 

In  this  section,  we  examine  the  effects  of  the  shift,  X5,  produced  by  inserting  an  optical 
window  at  an  angle  with  respect  to  die  optical  axis  after  the  acousto-optic  cell.  The  light 


distribution,  immediately  after  the  acousto-optic  cell,  is  represented  by 


(4.20) 


where  m  is  the  modulation  index,  a(x)  is  the  aperture  function,  and  O},  equals  where  61,  is  the 


illumination  angle  on  the  acousto-optic  cell.  This  light  distribution  passes  throu^  the  optical 
window,  which  translates  it  in  the  direction  of  propagation,  or  the  positive  x-direction,  of  the 
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acousto-opdc  cell  The  shifted  signal  distributioD  is  given  by 

f  f  *\  ■  f  \(  ^  j‘^^cU-^2  ^  >2***<*-'5) 

/+(x,f)s;ma(x-JC5)j  f---i - v  2  v 

\  2  V  )  (4.21) 

We  nowmni  our  attention  to  the  reference  beam,  which  is  a  function  of  x  only.  The  lens 

creates  an  effective  point  source  a  distance  Xq  from  the  optical  axis,  die  reference  beam  is  therefore 

represoited  by 

rix). 


.  TjC-Xq 
i  =  smc  — , 

.  ^  . 


(4.22) 


where  do  is  die  size  of  the  point  source.  The  reference  beam  and  die  signal  beam  combine  at  the 
beam  combiner,  and  are  imaged  onto  two  wideband  photodetectors. 

Since  the  photodetectors  are  in  an  image  plane  of  the  acousto-optic  cell,  the  signal  beam  and 
die  reference  beam  are  represented  by  reversed  coordinates.  We  let  u  be  the  coordinate  at  the 
photodetectors.  The  signal  light  distribution  at  the  photodetectors  is  represented  by 


Siimlarly,  the  reference  beam  at  the  photodetectors  is  represented  by 

.  .  .  , 
r(«)  s  smcf 

The  intensity  of  the  combined  beams  is 


(4.23) 


-U-Xq) 


(4.24) 

(4.25) 


l(u,0  =  +|r(w)j^  +2Re{4(«,/)r*(u)}, 

where  we  are  interested  in  only  the  cross-product  term,  since  the  others  are  removed  by  a  banc^ass 

filter.  The  photodetector  integrates  the  light  intensity  over  the  output  plane  such  that  g(t)  is  given 
by 

where  P(u)  is  the  photodetector  aperture  function.  In  this  analysis,  we  assume  that  P(u)  equals 
over  the  surface  of  the  photodetector  The  third  element  is  the  cross-produa  tenn,  which  is  given 
by 

g3(0  =  2Re|r  f^iu,t)r* iu)diA. 

t*'—  J  (4.27) 
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We  substitute  equations  (4.23)  and  (4.24)  into  (4.27)  to  get 


g3(/)a2Re 


(4.28) 


We  take  the  temporal  terms  outside  the  integral  and  sinq>lify  §3  (t)  to  get 

(4.29) 

-  We  tqjproximate  the  sine  function  as  a  delta  function,  as  we  did  in  Section  3.2.2,  to  get 

/  \  ^T>  /  •./  7*  (-U-X5)'\  ^  ~  gf 

^3(0  =  2Rej /m«  '>  0:5)1!  r-— - ^  S{-u- x^pu  \ 


(4.30) 

We  invoke  the  sifting  property  of  the  delta  function  to  get 

^3 (r)  =  ma{-XQ  - 0:5 )s^r  - ^ + jcos^2;?rc^r - y  +  j j- 

(4.31) 

where  we  realize  diat  the  ouqiut  of  the  photodetector  is  proportional  to  the  input  signal,  as  in 
Section  3.2.2,  but  with  an  added  time  delay  of  X5/V.  We  note  that  the  carrier  is  also  shifted  by  a 
phase  factor  which  is  directly  proportional  to  the  amount  of  the  physical  shift,  but  this  does  not 
affect  the  recovery  of  the  signal.  Finally,  a  scaling  factor  is  introduced  because  the  aperture 
function  is  probed  at  a  slightly  different  place  then  when  the  light  distribution  is  un.shifted, 
however,  the  difference  in  amplitude  is  minute  and  does  not  significantly  affect  the  current  output 
4.3  Measurements  on  the  optical  window 

Before  inserting  the  optical  window  into  the  Presort  Processor,  we  evaluated  the  optical 
window’s  baseline  performance.  An  830  nm  laser  produces  10  mW  of  power  incident  on  the 
optical  window.  The  transmitted  power  is  94.6  percent  or  9.46  mW  of  the  incident  light  The 
reflected  power  off  the  first  surface  (uncoated)  is  390  ^W,  and  the  power  reflected  off  the  second 
surface  (coated  widi  Hebbar  1016)  is  16.6  ^W.  We  calculate  the  percent  of  reflected  li^t  off  the 
first  surface  as 


0.390mW 


100  =  3.9%. 


lOmW 

Since  we  showed  in  Section  4.2.3  that  the  theoretical  power  reflected  off  the  first  surface  is 
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A  percent,  the  measured  value  and  the  calculated  value  are  in  close  agreement  The  power  reflected 
off  the  second  surface  is  0. 166  percent  of  the  incident  power,  which  compares  to  a  specification 
(from  Melles  Griot)  value  of  0.37  percent  If  we  add  up  the  percentages  calculated  above,  die  total 
l)ower  reflected  and  transmitted  is  98.77  percent  We  attribute  die  1.23  percent  of  light 
unaccounted  for  as  the  amount  the  optical  coating  absorbs. 

Returning  to  Section  4.2.3,  we  compare  the  measured  value  of  the  fringe  pattern  to  the 
theoretical  value.  The  amount  of  light  reflected  off  the  first  and  second  surfaces  are  3.9  and  0.166 
percent  of  die  incident  light,  respectively.  We  calculate  the  first  surface  Fresnel  reflection 
coefficient  as 

rpi=Va039  =0.1975. 

Similarly,  we  confute  the  second  surface  reflection  coefficient  as 

rp2=V0.00166  =  0.0407. 

We  calculate  b  as 


rp2  0.0407  « 

ft a  T— :::  =  0.206L 

Fpi  0.01975 

Hence,  we  now  calculate  the  inwnsity  in  the  Fourier  plane  as 

1(a)  -l+b^+ 2bcos(27tca^X 

or,  after  substituting  the  value  of  b,  we  get 

1(a)  =  1  +  0.0425  +  0.4122cos(2;rocc3). 

’  Therefore  a  30  dB  narrowband  interferer  will  range  between  28  dB  to  31.64  dB  depending  on  die 
value  of  a;  the  difference  between  the  experimental  results  and  die  calculated  values  is  attributed  to 
the  anti-reflection  coating  value  being  a  worst  case  specification.  Hence,  the  amount  of  distortion 
is  small  and  therefore  tolerable  for  our  application. 

'4.4  Insertion  of  the  optical  window  into  the  Presort  Processor 

Figure  4.7  shows  the  amplitude  and  phase  response  of  the  Presort  Processor  before  the 
optical  window  is  inserted.  The  amplitude  response  is  flat  to  within  3  dB,  except  for  die  six 
notches,  denoted  by  die  numbers  one  through  six,  in  the  passband.  These  notches  are  caused  by 
six  disfunctioning  spatial  li^t  modulator  elements;  the  elements,  which  have  a  one-to-one 
correspondence  to  frequencies  in  MHz,  are  266, 308, 348, 355, 421,  and  446.  The  bottom  part  of 
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the  trace  shows  the  phase  response  of  the  Presort  Processor  Notice  that  there  is  45  degrees  of 
phase  variation  over  the  200  MHz  bandwidth. 

We  dien  inserted  die  optical  window  into  the  interferometer  of  the  Presort  Processor  at 
approxiinately  a  twenty  degree  angle  and  the  anplitude  response  remained  unchanged.  Figure  4.8a 
shows  that  the  phase  response  of  the  Presort  Processor  is  linear,  note  the  sawtoodi  phase 
waveform  is  caused  by  die  network  analyzer  having  a  modulo  2n  scale.  The  linear  phase,  which 
corresponds  to  a  time  delay,  is  caused  by  the  0.558  mm  shift  of  the  optical  beam  in  die  direction  of 
propagation  of  the  acousto-optic  cell. 

We  null  out  the  linear  phase  by  introducing  an  electrical  time  delay  of  0.1328  |is;  this  time 
delay  corresponds  to  a  physical  shift  in  the  beam  of  0.558  mm.  Althou^  the  theoretical  value 
calculated  in  Section  4.3.2,  based  on  a  nominal  window  angle  of  20  degrees,  is  0.684  mm,  the 
experimental  result  given  here  implies  that  the  window  angle  was  actually  inserted  at  an  angle  of  17 
degrees. 

Figure  4.8b  shows  die  phase  of  the  Presort  Processor  after  the  optical  window  is  inserted 
'-and  after  the  linear  phase  is  nulled.  Notice  that  die  phase,  when  compared  to  the  phase  response  in 
Figure  4.7,  is  quadratic.  The  change  in  the  phase  response  is  attributed  to  the  path  difference 
obtained  when  the  optical  window  is  placed  at  an  angle  in  the  interferometer  of  die  Presort 
Processor.  Recall  that  die  focus  of  the  lens  in  the  upper  branch  should  be  an  equal  distance  from 
the  beam  combiner  as  the  signal  acousto-optic  cell.  However,  when  the  optical  window  is  placed 
inside  the  interferometer  at  an  angle,  die  padi  length  in  the  lower  branch  increases  by  die  cosecant 
of  the  angle  die  ray  inside  the  material  makes  with  respect  to  the  normaL  Therefore,  the  focal  plane 
of  the  lens  in  the  upper  branch  of  die  interferometer  and  the  acousto-optic  cell  are  no  longer  the 
same  distance  from  the  beam  combiner  The  lens  in  the  upper  branch  can  be  moved  closer  to  the 
first  beamsplitter,  to  correct  the  path  difference,  thus  increasing  the  distance  from  the  focus  to  the 
beamcombinec  Other  solutions  are  to  place  an  optical  window,  with  the  same  thickness  and  at  the 
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same  angle,  in  the  upper-branch,  or  to  move  the  acousto-optic  cell  closer  to  the  beam  combinet 

Figure  4.9a  shows  an  expanded  view  of  the  an^litude  re^onse  before  the  optical  window 
is  inserted  in  difrinterferometer  of  the  Presort  Processor:  Notice  that  the  notch  is  a  fraction  of  a 
MHz  to  die  right  of  330  MHz.  Figure  4.9b  shows  the  anqilitude  response  when  the  optical 
window  is  inserted.  The  notch  has  moved  to  the  tight  of  the  vertical  line  by  t^proximately  one- 
sixth  of  a  megahertz;  the  shift  is  caused  by  the  wedge  on  the  optical  window.  The  physical  shift  of 
.  die  spectrum  does  not  affect  the  performance  of  the  Presort  Processor  because  die  spatial  light 
modulators  are  moved  a  corresponding  amount,  using  their  micrometer  adjustments,  to 
accommodate  the  shifted  spectrum. 

Figure  4.10  shows  the  phase  response  when  a  second  optical  window  is  positioned  in  the 
upper  branch  of  the  interferometer,  at  the  same  angle  as  die  first  optical  window.  Figure  4. 10 
shows  that  the  phase  response  is  similar  to  the  phase  response  before  the  optical  windows  are 
inserted,  as  shown  in  Figure  7,  therefore  correcting  the  path  length  mismatch  between  the  reference 
beam  and  the  signal  beam.  The  sudden  jumps  in  phase  are  due  to  the  faa  that  six  of  the  spatial 
light  modulator  elements  are  not  functioning,  causing  six  notches  in  the  spectrum.  Since  the  noise 
floor  is  greater  than  the  notched  depth,  the  random  phase  jumps  occur  at  die  disfunctioning  spatial 
light  modulator  element  locations. 

As  a  completion  to  this  exercise,  we  modulated  a  pulse,  passed  it  through  the  Presort 
Processor,  demodulated  the  pulse,  and  displayed  it  on  an  oscilloscope.  We  performed  this 
experiment  before  the  primary  optical  window  is  inserted,  after  the  primary  optical  window  is 
inserted,  and  after  the  corrective  optical  window  is  inserted.  Figure  4. 1 1  shows  the  results  of  this 
experimenL  The  top  trace  is  without  the  primary  optical  window,  the  second  is  with  the  primary 
optical  window,  and  the  diird  is  after  die  corrective  optical  window  is  inserted. 

Note  that  the  rise-time  and  fall-time  of  the  pulse  increases,  hence  spreading  the  pulse,  from 
the  top  trace  to  the  middle  trace.  The  spreading  of  the  pulse  is  readily  understood  by  considering 
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the  path  difference  created  by  the  optical  window.  Since  the  path  length  from  the  acousto>optic  cell 
to  die  beam  combiner  is  greater  than  the  padi  length  from  the  focus  to  the  beam  combiner,  die  focal 
plane  of  the  lenain  die  upper  branch  is  effectively  outside  the  plane  of  the  acousto-optic  cell  The 
light  goes  through  its  focus  and  to  a  Fresnel  plane,  the  function  at  the  Frcsnel  plane  can  be  thought 
of  as  a  larger  point  source,  but  not  a  delta  function,  that  probes  the  signal  in  the  acousto-optic  cell. 
We  approximate  the  larger  source  as  a  rect  function  so  diat,  at  the  photodetector  plane,  the 
convolution  of  the  signal  with  the  rect  function  causes  the  pulse  to  spread,  as  second  trace  shows. 
From  the  third  trace,  we  notice  that  when  the  second  optical  window  is  inserted  in  the  upper  branch 
of  the  interferometer  to  bring  the  reference  focal  plane  into  coincidence  with  the  acousto-optic  cell, 
the  rise-time  and  fall-time  return  to  normal,  and  the  passband  of  the  pulse  is  restored  to  its  original 
value. 

4.5  Calculation  of  the  Fourier  lens  focal  length 

Rrst  we  calculate  the  length  of  the  D  scries  linear,  1024  element  CCD  photodetector  array 

that  is  to  be  illuminated.  We  choose  the  Fourier  transform  to  illuminate  800  photodetector 

elements,  which  are  13  micrometers  wide,  so  that  a  frequency  resolution  element  is  sampled  at  a 

four-to-one  ratio.  Hence,  the  array  length  illuminated  is 

A^=800(13^)=10.4  mm. 

The  spatial  variable  in  the  Fourier  plane  is  related  to  the  parameters  of  the  signal  and  the 
acousto-optic  cell  by 


V  (4.32) 

where  Af  is  the  bandwidth  of  the  signal,  X  is  the  wavelength  of  the  light,  v  is  die  acoustic  velocity 
within  the  cell,  and  F  is  the  focal  length  of  the  Fourier  lens.  We  have  Af=200  MHz,  X=830  nm, 
vs4200m/s,  and  10.4  mm  in  this  design.  We  reanange  equation  (4.32)  to  calculate  the 
required  focal  length  of  the  lens  as  F=263.1  mm. 


The  beam  height  at  the  acousto-optic  cell  is  250  micrometers  and  the  height  of  the  elements 
on  the  photodetector  array  is  13  micrometers.  Therefore,  to  collect  all  the  li^t,  a  demagnification 
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on  the  order  of  20  must  be  accomplished  by  using  a  cylindrical  lens. 

4.5.1  Fourier  lens  calculation 

We  first  calculate  the  focal  lengths  of  the  lenses  used  for  the  Fourier  transform,  and  we  then 

calculate  the  focal  length  of  die  imaging  lens.  The  closest  focal  length  lens  available  ofif-che-self  is 

260  mm.  To  accurately  illuminate  the  photodetector  array,  we  must  therefore  use  two  lenses. 

The  equivalent  power  of  two  lenses  is 

^eq  =  + 1^2  -  duKiKi,  (4.33) 

where  Kj  and  K2  are  the  powers  of  die  two  lenses,  and  d|2  is  the  distance  between  the  lenses;  the 

power  of  a  lens  is  the  reciprocal  of  its  focal  length.  We  use  two  available  lenses  whose  focal 

lengths  are  Fi=350  mm  and  F2=750  mm.  To  achieve  an  equivalent  focal  length  of  263.1  mm,  we 

calculate  die  distance  between  the  lenses  as  d]2=101.9  mm. 

Figure  4.12  shows  the  top  view  of  the  optical  system.  Using  ray  tracing  techniques,  h3, 

the  output  height,  is  derived  in  terms  of  the  input  angle,  the  distances  between  the  lenses,  and  the 

power  of  the  lenses.  The  height,  hj ,  at  the  first  lens  is 

*1  =  ^1^1’  (4.34) 

where  d^  is  the  distance  from  the  input  plane  to  the  lens,  and  6^  is  the  cutoff  angle.  We  calculate 
the  angle,  62,  shown  in  Figure  4.12,  using  the  refraction  equation 

^2  =  -  Ai^'i,  (4  35) 

where  is  the  power  of  the  first  lens.  We  continue  with  the  ray  tracing  technique  to  get  the 


(4.36) 

(4.37) 

(4.38) 


following  equadons: 

h2  —  h^  +  02^29 
^  =  ©2  “^2^2* 

^  ^  +  03^5, 

where  the  various  parameters  are  shown  in  Figure  4.12.  We  substitute  equations  (4.34)  through 
(4.37)  into  (4.38)  to  get 

^  =  &i(di  +  d2+di -didiKx -did^Kj -did^K^^,  (4.39) 

Since  d2,  6  j ,  the  height,  h3,  in  the  output  plane,  and  the  power  of  the  lenses  are  known,  we  can 


express  ds  as  a  function  di : 
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’■  A(l-<kfC^-d,K^)  ■ 

Since  we  have  the  output  distance  as  a  function  of  input  distance,  we  can  design  the  optical  system 
based  on  total  length  considerations. 

4.5.2  Design  and  placement  of  the  Fourier  transform  lenses 

Figure  4.1  shows  a  top  view  of  the  interferometer  of  the  Preson  Processor  The  distance 
from  the  acousto-optic  cell  to  the  optical  window  is  105  mm.  We  choose  tiie  distance  from  the 
acousto-optic  cell  to  the  first  lens  as  150  mm,  v^ch  places  the  lens  45  millimeteTS  outside  of  the 
interferometer.  The  lenses  in  the  optical  design  are  placed  outside  of  the  interferometer  so  they  do 
not  disturb  the  operation  of  tiie  Presort  Processor. 

Using  the  relationship  developed  in  the  last  section  between  the  input  and  output  distance, 
we  calculate  the  distance,  d3,  from  tiie  second  Fouria-  transforming  lens  to  the  photodetector  array. 
We  have  that  h3=5.2  mm,  0^=0.01976  radians,  di*150  mm,  d2*101.9  mm,  Fi=350  mm, 
F2=750  mm,  and  Fgq=263.1  mm.  We  substitute  the  parameters  into  equation  (4.40)  to  get  that 
d5=187  mm.  Therefore  the  total  length  of  the  optical  system  is  di+d2-Kl3=439  mm. 

4.5.3  Image  direction  calculations 

Figure  4.13  shows  the  image  direction  of  the  optical  system.  Since  the  output  spot  size, 
the  input  spot  size,  the  distances  of  the  Fourier  transform  lenses,  and  the  powers  of  the  first  two 
lenses  ate  known,  the  unknown  parameters  in  the  system  are  the  distances  d3,  d4,  and  F3.  We 
derive  d3  as  a  function  of  d4,  using  ray  tracing  techniques,  to  get 


_  —0444  ~  ^2  ~  d^d2^i ) 


(4.41) 


\~r.-rx; 

where  ail  the  parameters  except  for  d3  and  d4  are  known.  We  substitute  in  tiie  values  of  the 
parameters  to  get 


43  =  65.4144  -  708.24. 

From  tiie  Fourier  transform  design  we  know  that 

43  +  44  =187mni. 


(4.42) 


(4.43) 
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We  solve  equation  (4.43)  for  and  substitute  its  value  into  equation  (4.42)  to  get 

-<^4  + 1 87mni  *=  65.41df4  -  708.24.  (4.44) 

We  solve  equation  (4.44)  to  get  Aat  d4=13.S  mm.  We  now  substitute  the  value  of  d4  into  equation 

(4.43)  to  get  that  d3=:173.5  mm.  From  the  transfer  equation  we  have  that 

^  ^4  _  ^2  ~  ^^2  ~  ^4 

hi  -^4^4  ’  (4  45) 

where  all  the  parameters  in  the  equation  are  known.  We  find  that  K3s0.075317  mm,  so  diat  the 
focal  length  of  F3  is  1327  mm.  Melles  Griot  has  a  cylindiical  lens  with  a  focal  length  of  12.7  mm 
in  stock,  which  gives  sufficient  light  on  the  photodetector  array  for  die  integration  time  we  are 
using.  Summing  up,  we  have  F3=12.7  mm,  d3=173.5  mm,  and  d4=13.5  mm. 

4.6  Introduction  to  the  post-processing  electronics 

The  major  purpose  of  this  study  is  to  show  that  the  switching  of  the  spatial  light  modulator 
elements  can  be  automated  in  an  open-loop  fashion.  To  achieve  this  goal  we  implemented  an 
optical  power  spectrum  analyzer  to  read  out  die  power  spectrum  into  post-pnxressing  electronics; 
the  post-processing  electronics  were  constructed  using  wire-wrapping  techniques.  The  post¬ 
processing  electronics  locate  the  narrowband  inteiferers  in  the  spectrum,  and  output  a  location 
corresponding  to  the  spatial  light  modulator  element  where  the  narrowband  interferer  occurs.  A 
thresholding  scheme  is  used  to  achieve  this  operation,  along  with  a  feedback  loop,  to  adaptively  set 
the  threshold. 

4.6.1  Description  of  the  post-processing  electronics 

Figure  4.14  shows  a  block  diagram  of  the  post-processing  electronics.  The  power 
spectrum  of  the  signal  is  detected  by  the  CCD  photodetector  array,  and  is  sequentially  clocked  into 
the  anaiog-to-digital  converter,  which  converts  the  spectrum  to  a  digital  number  stream  The  digital 
numbers  are  compared  to  a  threshold  which  is  adaptively  set  by  the  post-processing  electronics,  or, 
alternately,  the  threshold  can  be  set  by  the  operatar  If  the  digital  number,  which  corresponds  to 
the  magnitude  of  the  spectrum  is  greater  than  the  tiireshold,  we  decide  tiiat  a  narrowband  interferer 
is  present;  conversely,  if  the  digital  number  is  less  than  the  threshold,  we  decide  that  no  interferer  is 
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present  The  location  of  die  pixels  in  the  photodetector  array  whidi  exceed  the  threshold  are 
counted  and  clocked  into  the  excising  routine.  The  excising  routine,  which  uses  a  look-up  table 
format,  outputs  die  spatial  light  modulator  element  location  information  required  to  perform  die 
switching. 

4.6.2  Adaptation  mode 

In  die  adaptation  mode  of  operation,  the  post-processing  electronics  dynamically  adi^ts  to  a 
steady-state  solution.  Sqipose  that  switches  one  and  two  of  the  block  diagram  shown  in 
Figure  4. 14  are  closed,  while  switch  three  is  open.  In  this  mode  of  operation,  the  operator  inputs 
the  number  of  narrowband  interferers  to  be  excised,  and  the  threshold  search  feedback  loop  sets  a 
threshold  level  where  the  number  of  narrowband  interferers  input  by  die  operator  are  excised.  We 
use  an  exan^le  to  illustrate  the  adaptation  mode  of  operation. 

Figure  4.15  shows  a  spectrum,  distorted  by  six  narrowband  interferers.  Suppose  that  the 
operator  sets  the  number  of  narrowband  interferers  to  be  excised  at  four.  The  first  time  the 
spectrum  is  read  out,  the  system  automatically  sets  the  threshold  to  the  largest  value.  The  threshold 
search  depends  on  whether  the  number  of  narrowt^ind  interferers  input  by  the  operator  exceeds,  is 
less  dian,  or  is  equal  to  the  number  of  narrowband  interferers  that  die  threshold  excises  fiom  the 
spectrum.  Since  none  of  the  narrowband  interfeim  exceed  the  threshold,  none  of  the  spatial  light 
modulator  elements  are  switched.  The  number  of  narrowband  interferers  exceeding  the  threshold, 
which  is  zero,  is  then  compared  to  the  number  of  interferers  the  operator  wants  to  excise.  Since 
this  number  is  less  than  the  number  set  by  the  operator,  the  search  routine  continues.  The  second 
time  the  array  is  clocked  out,  the  post-processing  electronics  sets  the  threshold  to  the  stqi  below 
the  initial  value,  where  diere  are  16  equally  spaced  steps  across  the  range  of  an  eight-bit  word.  The 
number  of  narrowband  interferers  exceeding  this  threshold  are  counted,  conqiared  to  the  number  of 
narrowband  interferers  the  operator  wants  to  excise  and  the  ^propriate  spatial  light  modulator 
elements  are  switched  to  remove  those  narrowband  interferers  that  exceed  the  threshold. 
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For  each  readout  of  the  anay,  the  number  of  naxrowband  interferers  exci^  is  compared  to 
.the  number  of  narrowband  interferers  the  operator  desires  excised,  and  sets  a  new  direshold 
depending  on  whether  the  number  of  narrowband  interferers  excised  is  greater  than,  or  less  than 
the  desired  number.  When  the  desired  number  of  narrowband  interferers  excised  is  met,  the 
threshold  is  stationary  until  a  new  value  of  narrowband  interferers  is  input  by  the  operator  or  until 
f  some  narrowband  interferers  disappear  or  new  ones  appear  in  the  spectrum. 

Operator  interaction  mode  and  visual  display 

In  the  operator  interaction  mode,  switches  one  and  two  are  open  while  switch  three  is 
closed.  The  spectrum  is  sequentially  clocked  off  the  array,  converted  to  digital  numbers  while  die 
operator  sets  a  threshold  to  which  the  values  of  the  spectrum  are  compared.  Figure  4. 15  shows  a 
spectrum  of  a  signal  that  is  distorted  by  six  narrowband  interferers.  The  operator  reduces  or 
increases  the  threshold  depending  upon  how  many  of  the  six  narrowband  interferers  he  wants  to 
excise.  If  he  prefers  to  have  only  four  excised,  the  position  shown  is  correct  Since  the  operator 
monitors  the  spectrum  and  the  threshold  via  the  display  in  this  mode  of  operation,  the  operator 
:  completes  the  feedback  loop. 

The  visual  display  enhances  the  operation  of  the  system  by  showing  the  power  spectrum  of 
the  signal  and  the  threshold  level  of  the  post-processing  electronics.  Figure  4.16  shows  a 
spectrum  on  the  oscilloscope  display  with  an  associated  threshold.  The  spectrum  of  the  signal  is 
tapped  directly  off  the  CCD  photodetector  array  and  feeds  channel  one  of  the  oscilloscope;  the 
display  is  triggered  from  chaimel  one.  The  threshold  feeds  chaimel  two  of  the  oscilloscope  which 
;  appears  as  a  horizontal  line  across  the  display.  Every  part  of  the  spectrum  appearing  above  the 
horizontal  line  is  excised  while  the  remaining  spectrum  is  undisturbed,  the  lone  interferer  in 
Figure  4.18  is  therefore  excised. 

4.6.4  Excising  algorithm 

The  excising  algorithm  is  a  versatile  pan  of  the  post-processing  electronics.  The  operator 


49 


decides  the  expected  type  of  spectral  environment,  and  programs  an  algorithm  tlut  models  the 
environment  into  the  post-processing  electronics.  The  excising  algorithm  uses  the  width  of  the 
interfeier,  which^s  counted  by  a  threshold  counter,  together  with  the  location  of  the  interferer 
relative  to  the  spatial  light  modulator  elements  to  decide  which  spatial  li^t  modulator  elements  to 
switch.  As  an  exan^le,  suppose  an  interferer  is  located  in  the  center  of  the  fourth  spatial  light 
modulator  element  and  has  a  width  of  one  spatial  light  modulator  element  If  the  operator  wants  to 
switch  only  the  spatial  light  modulator  element  on  which  die  interferer  occurs,  the  excising 
algorithm  ou^uts  the  address  which  corresponds  to  spatial  li^t  modulator  element  number  four 
However,  if  die  operator  desires  to  switch  the  neighboring  spatial  light  modulator  elements,  he  can 
change  the  algorithm  by  reprogramming  the  post-processing  electronics. 

The  CG^  photodetector  array  samples  the  spatial  light  modulator  elements  at  a  ratio  of  four 
CCD  photodetector  elements  per  one  fin^uency  resolution  element,  equal  to  a  spatial  light 
modulator  element  width,  so  that  the  electronics  can  decide  where  the  interferer  is  located  relative  to 
the  spatial  light  modulator  elements;  the  excising  algorithm  is  stored  in  electrically  programmable 
read-only-memories  that  are  ultraviolet  erasable.  The  excising  algorithm,  used  for  our 
experiments,  searches  for  interferers  that  have  a  width  of  one  frequency  resolution  element  If  the 
interferer  is  located  directly  on  a  spatial  li^t  modulator  element,  that  element  is  switched  plus  the 
neighboring  elements;  no  more  than  three  elements  are  ever  switched  for  one  interferer.  If  the 
interferer  is  principally  located  on  two  spatial  light  modulator  elements,  then  both  ^atial  light 
modulator  elements  are  switched,  plus  the  spatial  light  modulator  element  that  is  the  nearest 
neighbor  to  the  element  on  which  the  most  energy  of  the  interferer  is  incident  If  the  interferer  is 
located  directly  between  two  spatial  light  modulator  elements,  they  are  switched,  plus  the  spatial 
li^t  modulator  element  that  is  the  nearest  neighbor  to  the  lower  number  spatial  light  modulator 
element  upon  which  the  interferer  is  incident  The  two  primary  spatial  light  nxxlulator  elements 
plus  a  nearest  neighbor  needed  to  be  switched;  the  nearest  neighbor  could  be  tiie  hi^er  number  or 
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die  lower  niunber  spatial  li^  modulator  elemenL  To  verify  the  operation  of  the  post-processing 
electronics  a  variety  of  interferer  scenario’s  were  implemented  and  the  results  are  given  in  the 
Section  4.6. 

4.7  Experimental  results 

Our  primary  design  goal  is  to  sort  the  spectrum  of  the  signal,  locate  the  narrowband 
interferers,  and  supply  the  spatial  light  modulator  element  addresses  in  less  than  10  milliseconds. 
.The  post-processing  electronics,  which  were  wire-wr^ped  and  delivered  to  RADQ  sorts  the 
entire  spectrum,  including  the  224  CCD  elements  that  are  not  illuminated,  and  outputs  the  spatial 
light  modulator  addresses  for  a  particular  scan  time  in  5.4  milliseconds;  we  tl^fore  have  met  our 
primary  design  goal.  Figure  4. 17  shows  the  ejq)erimental  setup  used  to  test  the  electronics.  The 
•  laser  beam  is  e^anded  by  lenses  and  and  passes  through  a  polarizer  which  is  used  to 
control  the  beam  intensity.  The  beam  is  incident  u{x>n  a  phase-grating  which  produces  several 
diffracted  orders.  The  Fourier  transform  of  the  phase  grating  is  formed  by  lens  L3  at  the  plane  of 
the  CCD  photodctcctor  array,  resulting  is  numerous  spots,  which  simulate  narrowband  interferers, 
along  the  length  of  the  array.  The  addresses  from  the  electronics  are  decoded  and  fed  to  a  10 
element  LED  display,  which  are  used  to  observe  which  spatial  light  modulators  are  switched. 

We  test  the  adaptation  mode  of  operation  by  inputting  the  number  of  interferers  to  be 
excised,  and  showing  that  the  post-processing  electronics  adapt  to  meet  out  specification.  The  first 
test  consists  of  one  interferer  dominating  the  spectrum  and  we  vary  the  number  of  interferers  to  be 
excised.  We  input  zero  interferers  causing  the  threshold  to  remain  at  the  maximum  value.  We 
•increase  the  interferers  desired  excised  to  one  and  the  threshold  decreased  to  the  steady-state  value, 
as  Figure  4.18  shows,  therefore  ouqmtting  the  ccnrect  addresses.  Since,  the  location  of  the 
interferer  is  known,  we  confirmed  that  the  correct  spatial  light  mndulatnr  elements  are  switched. 

The  second  experiment  is  to  have  two  interferers  of  different  strengths  while  the  number  of 
interferers  desired  excised  is  varied.  If  we  require  that  no  interferer  be  excised,  the  threshold 
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remains  at  die  maximum  value.  Figure  4.19  shows  the  result  when  we  require  that  one  interferer 
be  excised;  the  electronics  successfully  finds  the  coirect  steady-state  threshold.  When  we  require 
that  two  interferers  be  excised,  the  threshold  decreased  until  botii  inteiferers  were  successfully 
notched. 

In  dynamic  environments,  nairowband  interferers  are  not  necessarily  stationary  in 
anq}litude  or  frequency,  but  are  more  likely  to  be  agile  interferers  that  hop  through  the  spectrum 
according  to  a  pseudo-random  sequence.  To  simulate  an  intearfeier  hopping  through  the  spectnon, 
we  varied  tiie  location  of  the  interferers  on  the  CCD  array,  and  monitored  the  ^atial  light 
modulator  elements  that  were  switched.  We  changed  the  position  of  tiie  interferer  by  sliding  the 
CCD  photodetector  array  horizontally  while  tiie  electronics  processed  the  new  data.  The 
electronics  tracked  tiie  interferer  to  its  new  location  and  ouqmt  the  addresses  to  switch  the  correct 
spatial  light  modulator  elements.  We  noticed  that  it  took  the  post-processing  elecminics  at  most 
two  scan  times  to  successfully  provide  the  address  needed  to  notch  the  interferex: 

Figure  4.20  shows  the  final  interferer  scenario  used  to  test  the  electronics.  There  are  four 
total  interferers  with  two  having  equal  aiiqilitudes  and  the  remaining  two  having  different 
strengths.  The  purpose  of  this  experiment  is  to  show  that  a  threshold  can  be  set  where  we  might 
only  want  to  excise  the  three  largest  interferers.  We  manually  varied  the  threshold  until  the  position 
shown  in  Figure  4.20  was  reached.  The  electronics  ouqiut  the  correct  addresses  to  successfully 
notch  the  three  interferers,  while  leaving  the  smallest  interferer  untouched.  We  performed  the  same 
battery  of  tests  as  performed  in  the  adaptation  mode  with  the  electronics  in  the  operator  interaction 
mode.  In  all  cases,  the  addresses  to  the  correct  spatial  light  modulator  elements  were  output  We 
have  therefore  con^leted  the  testing  of  the  electronics,  and  now  nx)ve  ahead  to  the  experimental 
results  obtained  on  the  Presort  Processor 
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5.  Pulse  experiments 

In  diis  ch^ter  we  explore  the  effect  of  removing  frequency  componoits  of  a  pulse  in  the 
Presort  Processor.  We  describe  the  results  of  coiiq)uter  simulations  and  e}q)eriments  in  which  we 
remove  spectral  components  of  a  pulse  at  an  intermediate  frequency,  translate  it  to  baseband,  and 
characterize  the  distortion;  these  experiments  utilize  the  multi-channel  acousto-optic  cells  in  the 
Presort  Processor  to  perform  the  frequency  removal.  We  study  how  this  frequency  removal  relates 
to  the  final  pulse  shape  at  the  photodetector,  using  the  root  mean  square  error  as  a  measure  of  the 
distortion.  Finally,  we  add  a  narrowband  interferer  to  the  pulse,  excise  it,  and  characterize  the 
distortion  after  die  pulse  is  translated  back  to  baseband. 

5.1  Signal  excision 

In  the  early  1980’s  Erickson  [1-3]  perfoniMd  analyses  and  experiments  dealing  with  signal 
excision.  He  used  a  Mach-Zehnder  interferometer  with  an  acousto-opdc  cell  in  one  branch  and 
beam  shaping  optics  in  the  other  branch  to  perform  his  experiments.  His  main  interest  was  in 
finding  the  notch  shape  and  depth  in  the  frequency  plane,  using  different  reference  beam  shapes 
and  a  variety  of  illuminations  on  die  acousto-optic  cell 

Bandstetter  and  Grieve  [7-8]  built  an  optical  notching  filter  that  uses  a  recursive  technique 
to  obtain  the  final  notch  depth.  They  use  a  programmable  spatial  filter  in  which  liquid  crystal 
fingers  block  the  undesired  light  Their  experiments  involves  passing  an  RF  pulse  through  their 
optical  excisor,  and  examiiiing  the  effect  their  system  has  on  the  pulse.  Narrowband  interferers  are 
then  added  to  the  RF  pulse,  excised  using  the  recursive  architecture,  and  dien  examined  at  the 
output  at  RF. 

5J  Computer  simulation  parameters  and  the  experimental  setup 

In  the  computer  simulation  model  that  simulates  the  I^esort  Processor,  we  use  a  Gaussian 
distribution  that  is  truncated  at  the  l/e^  points  to  illuminate  the  acousto-optic  cell  The  reference 
beam  is  represented  as  a  Gaussian  distribution,  also  truncated  at  the  l/e^  points,  at  the  frequency 
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plane.  We  let  the  contrast  ratio  of  the  ^atial  li^t  modulators  be  infinite,  and  assume  that  all 
signals  are  modulated  to  die  center  fiequency  of  the  signal  acoustoKiptic  celL 

In  die  first  set  of  expetiments  we  remove  certain  frequency  con^nents  of  the  spectrum  of 

a  pulse  at  an  intermediate  frequency,  and  analyze  the  resulting  distortion  on  die  pulse  after  it  is 

brought  to  a  baseband  fiiequency  range.  Figure  5.1  shows  a  block  diagram  of  the  experimental 

setup.  The  input  baseband  pulse  is  multiplied  by  a  cw  signal,  translating  the  spectrum  of  the  pulse 

to  the  intermediate  frequency  f^.  The  drive  signal  to  the  acousto-optic  cell  is  given  as 

f(t)  =  s(t)  cos(2ji  fgt). 

The  Presort  Processor  modulates  the  light  in  both  space  and  time,  and  Fourier  transforms  it  at  the 
fiequency  plane  containing  the  spatial  light  modulators.  Pan  of  die  spectrum  is  removed  by  the 
spatial  light  modulators  and  the  remaining  light  is  collected  onto  the  pbotodetectors.  The  output  of 
the  Preson  Processor  is  a  pulse  f^,  with  pan  of  its  spectrum  removed. 

The  signal  fixim  the  photodetectors  is  then  mixed  with  the  same  local  oscillator  at  fc, 
yielding  a  term  at  baseband  and  a  term  at  twice  the  intomediate  fiequency.  The  double  frequency 
term  is  rejected  by  a  low-pass  filter,  and  we  denote  the  remaining  distorted  pulse  by  d(t).  The 
distortion  is  measured  relative  to  the  output  signal,  b(t),  which  represents  the  signal  when  none  of 
the  spatial  light  modulators  are  switched. 

Figure  5.2  shows  a  mote  conqilete  block  diagram  of  the  experimental  setup.  The  HP 
8116A  pulse/function  generator  produces  a  pulse  which  is  fed  into  mixer  1.  An  intermediate 
firequency  is  generated  by  the  HP  8640B  signal  gmerator,  and  is  supplied  to  the  power  splitter 
The  number  2  output  of  the  power  splitter  is  connected  into  mixer  1  to  translate  the  baseband  pulse 
to  the  bandpass  of  the  Presort  Processor  The  signal  out  of  mixer  1  is  then  anq)lified  by  30  dB, 
using  a  one  watt  Amplifier  Research/Model  IWlOOO  amplifier.  The  amplified  signal  drives  the 
acousto-optic  cell  in  the  Presort  Processor 

The  signal  which  is  deteaed  at  the  photodetectors  is  connected  to  mixer  2.  Port  3  of  the 
power  splitter  feeds  the  local  oscillator  to  mixer  2.  The  ouq)ut  of  mixer  2  is  filtered  by  a  100  MHz 


low-pass  filter,  and  the  baseband  signal  is  displayed  on  a  Tektronix  1 1403  Digitizing  Oscilloscope. 
Measure  of  distortion 

There  are-vaiious  methods  for  measuring  signal  distortion.  In  die  fi:equency  domain,  filters 
are  characterized  by  parameters  such  as  passband  ripple,  overshoot,  and  undershoot  lime  domain 
distortion  if  often  measured  by  parameters  such  as  rise-time  and  fall-time.  For  our  experiments  we 
chose  to  use  the  root  mean  square  error  criterion  as  a  measure  of  distortioit  The  root  mean  square 
value  of  a  waveform,  v( rj,  is  given  by  _ 

where  T  is  the  duration  of  the  waveform.  From  Section  5.2  we  recall  that  the  undistoit^  signal  is 
denoted  as  b(t)  and  the  distorted  signal  is  given  as  d(t).  Therefore,  the  error  of  the  distorted  pulse 
relative  to  the  undistorted  pulse  is 

e(t)  =  b(t)-d(t). 

The  root  mean  square  error  is  given  by  _ 

or 

=  1^16(0- ./(Of*. 

In  the  following  experiments,  we  normalize  the  root  mean  square  error  to  the  root  mean  square 
value  of  the  undistorted  pulse  b(t). 

5.3.1  System  limitations  due  to  the  components 

Several  components  in  this  experimental  setup  limit  the  performance  of  the  system.  For 
example,  the  local  osdUator  has  a  tendency  to  jump  sixldenly  by  one-half  of  a  megahertz.  The 
pulse  generator  produces  pulse  shapes  that  have  large  overshoot  and  slow  fall  times.  The  mixers 
limit  die  operation  of  the  system  because  the  local  oscillator  input  has  a  low  maximum  input 
voltage.  If  the  maximum  input  voltage  level  is  exceeded,  the  mixer  electronics  are  overdriven, 
causing  spurs  which  are  similar  to  narrowband  interferers. 


Figure  5.3  shows  die  cleanest  signal  obtainable  fiom  the  experimental  setup  without 
overdriving  the  mixers;  this  pulse  is  what  a  receiver  would  see  on  a  one-shot  basis.  The  hei^t  of 
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the  pulse  is  7  millivolts  with  a  noise  level  of  2  millivolts,  which  conesponds  to  a  10.9  dB  signal- 
to-noise  ratio.  The  Presort  Processor  has  a  signal-to-noise  ratio  of  25  dB,  while  die  mixers  have  a 
signal-to-noise  ratio  that  is  much  less  than  25  dB,  resulting  in  the  mixers  limiting  die  overall 
system  perfcnnance.  Triggering  die  oscilloscope  to  display  die  wavefonn  was  dif&cult  due  to  the 
limitations  of  die  test  equipment  We  therefore  elected  to  overdrive  the  local  oscillator  input 
yielding  a  larger  signal  level  as  well  as  a  larger  noise  level  at  the  ouqiut  We  subsequently  averaged 
die  ouqiut  signal  256  times  to  reduce  the  noise  level  Since  the  averaging  is  equivalent  to  coherent 
detection,  the  signal-to-noise  ratio  attained  in  the  following  set  of  experiments  is  larger  than  die 
actual  signal*to-noise  ratio  at  the  ouqiut  of  the  system. 

5.3 J  The  removal  of  symmetric  frequency  components  on  a  78  ns  pulse 

In  dynamic  environments  there  are  numerous  emitters  that  may  render  a  communication 
system  useless.  We  want  to  excise  these  interferers,  using  frequency  plane  excision,  to  recover  the 
intended  message.  When  we  excise  the  narrowband  inierferers,  however,  part  of  the  energy  of  the 
intended  message  is  also  removed.  In  this  experiment,  we  study  the  effects  of  removing 
symmetric  frequency  components  of  a  78  ns  pulse  at  an  intermediate  frequency,  and  analyze  the 
output  pulse  shapes  at  baseband.  To  further  illustrate  these  experiments,  we  consider  the  sine 
function  in  Figure  5.4a  which  represents  the  FoutIct  transform  of  a  pulse.  Figure  5.4b  shows  that 
the  center  value  of  the  sine  function  is  removed  and  Figure  5.4c  shows  that  symmetric  frequency 
components  are  removed  ^ut  the  center  value.  The  removal  of  symmetric  frequency  components 
of  the  sine  function  is  equivalent  to  notching  a  sinusoid  whose  frequency  is  equal  to  the  frequency 
distance  in  Hertz  from  the  center  value  of  the  sine  function,  and  an  aii^litude  equal  to  the  height  of 
die  sine  function  at  that  particular  distance. 

Before  we  analyze  the  results  of  the  computer  simulation,  let  us  first  conjecture  what  we 
expect  to  obtaiiL  Since  the  frequency  components  are  removed  symmetrically,  we  should  expect 
the  error  to  be  largest  in  the  middle  of  the  band  and  roll  off  as  we  remove  frequency  conqionents 
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furdier  from  the  center  frequency  until  we  reach  die  first  null;  because  the  sine  function  has  90 
percent  of  its  energy  in  the  mainlobe  and  die  amplitude  of  the  sine  function  rolls  off  to  zero  as  we 
approach  die  fiistnulL  At  die  first  null,  the  error  should  be  a  minimum  and  then  oscillating  in 
magnitude  as  we  continue  notching  frequency  components  away  from  the  null. 

The  spectrum  of  the  pulse  in  the  conqiuter  simulations  is  centered  at  350  MHz  with 
sidelobes  as  described  in  the  previous  section.  Hrst,  we  simulate  the  switching  of  the  spatial  light 
modulator  element  350,  which  removes  the  central  value  of  the  sine  functiorL  Although  diree 
spatial  light  modulator  elements  are  switched  to  get  a  25  dB  notch  in  the  Presort  Processor,  the 
notch  positions  are  represented,  in  the  following  discussions,  by  the  frequency  of  just  the  central 
spatial  light  modulator  element  In  subsequent  experiments,  we  remove  frequency  components 
centered  at  348  and  352,  which  removes  symmetric  frequency  components  of  the  pulse.  The 
notches,  which  are  equal  to  three  spatial  li^t  modulator  elements  wide,  are  then  incremented  away 
from  the  center  of  the  spectrum  in  steps  of  one  spatial  light  modulator  element  'fherefore  the  next 
frequency  components  removed  are  centered  at  347  as  well  as  at  353. 

Figure  5.5  shows  the  computer  simulation  results  of  the  root  mean  square  error  as  a 
function  of  spatial  light  modulator  element  posidoa  The  root  mean  square  error  begins  at  0.50  and 
increases  to  0.68  at  352  MHz.  The  increase  in  error  is  due  to  the  removal  of  the  central  frequency 
components,  which  requires  that  only  three  spatial  light  modulator  elements  be  switched  off,  while 
the  removal  of  the  first  set  of  symmetric  frequency  components  requires  that  six  spatial  li^t 
modulator  elements  be  switched  off.  Since  the  symmetrical  notches  are  near  the  hipest  amplitude 
of  the  spectrum,  the  root  mean  square  error  increases  initially. 

Since  die  mainlobe  is  25.6  MHz  wide,  the  sine  function  reaches  a  minimum  at  362.8  MHz; 
we  therefore  expect  the  root  mean  square  error  to  be  a  minimum  ?t  362  or  363  MHz.  Figure  5.5 
shows  a  minimum  at  362  MHz  with  the  value  at  363  MHz  being  only  slightly  larger.  Therefore, 
the  computer  simulation  results  agree  with  our  expectations.  We  further  expect  the  toot  mean 
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square  error  to  follow  the  general  shape  of  a  sine  distribution,  as  shown  in  Figure  5.5,  since  we 
are  removing  die  energy  of  the  pulse  symmetrically  about  the  central  value. 

Figure  5.6  shows  the  experimental  results  of  symmetrically  removing  frequency 
coii^onents  of  the  spectrum  of  a  78  ns  pulse;  the  graph  shows  the  root  mean  square  error  as  a 
function  of  spatial  light  modulator  element  position,  which  has  a  one-to-one  mapping  to 
firequencies  in  MHz.  Note  that  the  pulse  is  centered  at  350  MHz  in  die  computer  simulations,  but 
in  our  experiments  the  pulse  is  centered  at  332  MHz.  Since  we  are  calculating  the  root  mean 
square  error  relative  to  the  undistorted  pulse,  as  symmetric  frequency  components  are  removed,  the 
results  are  the  same  since  we  remove  the  same  frequency  components  in  both  the  experimental  and 
computer  simulations.  Therefore,  to  make  the  comparison  between  experimental  and  the  computer 
simuladons  straightforward,  we  match  the  frequency  scale  of  the  experimental  results  to  those  of 
the  simulation.  The  root  mean  square  error,  at  the  center  frequency  of  350  MHz,  begins  at  0.29 
and  increases  to  0.37  at  352  MHz.  The  experimental  root  mean  square  error  is  smaller  than  the 
computer  simulation  results  because  infrnite  notch  depths  were  used  in  the  computer  simulations. 
We  expea  the  root  mean  square  error  to  be  minimized  at  362  or  363  MHz.  Since  Figure  5.6 
shows  a  minimum  at  362  MHz,  with  the  value  at  363  MHz  being  only  slightly  larger,  the 
experimental  results  agree  with  the  computer  simulation  results. 

We  expea  the  root  mean  square  error  to  increase  as  we  symmetrically  notch  through  the 
first  sidelobe.  This  phenomena  extends  to  367  MHz,  where  it  begins  to  decrease  instead  of 
increasing  as  we  expect  The  root  mean  square  error  does  not  increase  because  the  output  pulse 
displayed  on  the  oscilloscope  is  unstable.  The  instability  is  caused  by  the  tenqierature  dependence 
of  the  system  and  by  drifting  of  the  local  oscillator  frequency.  The  oscillator  frequency  was 
readjusted  at  367  MHz  to  obtain  a  good  input  pulse  shape,  and  the  renormalized  data  causes  the 
results  to  be  off  by  a  gain  faaot  With  this  in  mind,  we  recognize  the  root  mean  square  error 
would  trace  out  the  first  sidelobe  and  then  produce  a  null  at  spatial  light  modulator  element  374. 
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The  nx>t  mean  square  error  is  a  minimum  at  374  MHz,  and  increases  after  that  point  Another 
minimum  should  occur  at  386  MHz,  which  is  backed  tq>  experimentally  by  the  minimum  at 
386  MHz  in  Rgnre  5.6. 

Figure  5.7  shows  the  spectrum  of  die  experimental  pulse,  after  being  passed  through  the 
system,  but  before  being  demodulated  with  the  local  oscillator.  We  observe  that  the  mainlobe  and 
the  first  two  sidelobes  of  the  sine  function  are  strong  enough  to  recognize;  however,  we  observe 
that  the  third  sidelobe  is  close  to  the  noise  flooc  The  root  mean  square  error  after  386  MHz 
generally  increases,  widiout  following  the  sidelobe  pattern,  since  the  energy  level  of  the  signal  at 
die  output  is  close  to  the  noise  floor. 

Before  we  consider  the  experimentally  obtained  pulse  shapes  as  a  function  of  removing 
frequency  components,  we  predict  the  effect  on  the  pulse  shapes.  Fourier  series  analysis  shows 
that  all  signals  consists  of  weighted  frequency  harmonics.  The  frequency  harmonics  of  an  isolated 
pulse  consists  of  a  d.c.  value,  since  the  pulse  is  offset  from  die  horizontal  axis,  plus  a  continuum 
of  weighted  harmonics.  The  fundamental  harmonic  is  added  to  the  d.c.  value  to  get  a  pulse  diat  is 
rounded,  then  successive  harmonics  are  added  to  obtain  the  sharp  transitions  associated  with  the 
rise-time,  the  fall-time  and  the  flatness  on  the  top  of  the  pulse.  The  d.c.  value  and  the  fundamental 
reside  in  the  mainlobe  of  the  sine  function,  while  the  higher  hequency  components  reside  in  the 
sidelobes  of  the  sine  function.  We  dierefore  expect  that  the  removal  of  the  mainlobe  frequency 
components  will  produce  the  most  distortion  on  the  pulse,  and  the  removal  of  the  higher  frequency 
components  will  mosdy  affect  the  top  of  the  pulse,  die  rise-time,  and  the  fall-time. 

We  now  consider  the  shape  of  the  pulses  as  various  frequency  components  are  removed. 
Figure  5.8  shows  the  experimental  reference  pulse,  while  Figure  5.9  shows  the  experimental  pulse 
after  the  three  central  spadal  light  modulator  elements  are  switched.  Notice  that  a  dc.  value  has 
been  subtracted  frxim  the  pulse,  which  is  intuidvely  correct,  since  the  dc.  value  is  located  at 
350  MHz.  A  comparison  of  the  reference  pulse  with  the  distorted  pulse  reveals  that  the  siuqie  of 
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the  pulse  is  preserved,  but  with  a  d.c.  value  subtracted. 

We  now  examine  die  simulated  pulse  when  elements  344  and  356  are  switched  and 
compare  it  to  diecxperimental  pulse  when  die  same  symmetric  removal  of  frequency  components 
is  performed.  Hgures  5. 10a  and  5. 10b,  which  are  the  simulated  and  experimental  pulses, 
respectively,  show  that  part  of  the  fundamental  is  removed,  as  is  evident  by  the  sinking  of  the 
middle  of  the  pulse.  The  simulated  and  experimental  pulses  have  similar  sh^ies,  confirming  diat 
our  computer  model  is  correct 

When  we  remove  the  fiequency  components  in  die  first  sidelobe,  the  edges  of  die  pulse 
become  rounded,  as  Figures  5.1  la  and  5.1  lb  show,  because  we  are  removing  die  higher 
fi^equency  components.  The  fundamental  is  evident  by  the  hunp  in  die  center  of  the  pulse;  this 
simulated  pulse  is  with  spatial  light  modulator  elements  332  and  368  switched. 

A  similar  effect  happens  as  we  remove  fiequency  components  in  the  second  sidelobe  and 
beyond.  The  edges  on  the  pulse  become  rounded  and  ripple  is  evident  on  the  top  of  the  pulse  due 
to  the  removal  of  the  higher  frequency  components. 

5.3,3  The  removal  of  symmetric  frequency  components  on  a  468  ns  pulse 

The  efiects  of  excising  narrowband  intexferers  on  the  pulse  shape  is  a  function  of  the  pulse 
duration.  A  given  notch  width  represents  removing  a  large  fraction  of  die  signal  spectral 
components  when  the  pulse  is  wide,  because  the  spectrum  of  the  pulse  is  narrow.  To  illustrate  this 
effect,  we  repeated  the  frequency  removal  experiments  using  a  468  ns  pulse  as  the  ii^ut  signal. 
Figure  5.12  shows  the  spectrum  of  a  468  ns  pulse  afrer  it  is  mixed  to  baseband  by  the  local 
oscillator;  A  long  pulse  in  die  time  domain  produces  a  larger  signal-to-noise  ratio  in  the  frequency 
domain  than  a  short  pulse  because  the  energy  is  spread  over  a  smaller  region.  The  nulls  of  the  sine 
function  occur  every  2.14  MHz.  The  spectral  components  of  the  pulse  are  removed  in  the  same 
manner  as  for  the  78  ns  pulse,  but  in  this  case  there  are  only  two  spatial  li^t  modulator  elements 
per  sidelobe  and  that  the  mainlobe  spans  only  four  spatial  light  modulator  elements.  Since  90 
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percent  of  the  energy  in  a  sine  function  is  in  die  mainlobe,  die  root  mean  square  error  should 
decline  qmckly. 

Figure  S:13a  and  Figure  S.13b  show  die:  con^uter  simulation  and  experimental  results  of 
the  root  mean  square  error  as  a  function  of  spatial  light  modulator  element  position,  reflectively. 
The  root  mean  square  error  is  the  largest  when  we  switch  spatial  light  modulator  element  3S0.  The 
large  root  mean  square  error  is  due  to  the  concentration  of  energy  near  the  center  spatial  light 
'  modulator  element  By  switching  off  three  spatial  light  modulator  elements  we  are  eliminating  over 
50  percent  of  the  mainlobe.  We  note  that  the  shape  of  the  root  mean  square  error  of  the  computer 
simulations  are  similar  to  the  experimental  results. 

The  first  removal  of  symmetric  fiequcncy  components,  which  corresponds  to  elements  348 
and  352  being  switched,  is  less  traumatic  on  the  pulse  shape  as  confirmed  by  the  decrease  in  the 
root  mean  square  error.  The  decrease  in  the  root  mean  square  error  is  reasonable,  since  the  notches 
are  centered  two  spatial  light  modulator  elements  from  the  center  of  the  spectrum,  therefore  not 
seriously  afrecting  the  energy  in  the  mainlobe. 

When  we  move  the  symmetric  notches  one  more  increment  from  the  center  of  the  spectrum, 
we  notice  that  the  root  mean  square  error  drops  off  significantly  because  the  notches  now  are 
outside  the  mainlobe  of  the  sine  function.  The  root  mean  square  error  does  not  follow  the  shape  of 
a  sine  function  because  die  width  of  the  spatial  light  modulator  elements  are  large  with  respect  to 
the  distance  to  successive  nulls  on  the  sine  function. 

Figure  5.14a  and  Figure  5. 14b  show  the  computer  simulation  and  experimental  results, 
respectively,  of  the  pulse  when  we  switch  the  three  most  centered  spatial  light  modulator  elements. 
Notice  that  we  notched  the  d.c.  component,  plus  a  significant  amount  of  the  fundamental,  resulting 
in  a  differentiation  to  occur  on  the  pulse.  Also  note  the  similarity  in  the  simulation  results  when 
compared  to  the  experimental  results;  the  small  discrepancy  of  the  two  plots  is  caused  the 
removal  of  more  than  the  mainlobe  of  the  sine  function  in  the  experimental  results.  The  remaining 
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pulse  shapes  due  to  the  excisions  follow  a  pattern  similar  to  those  of  the  78  ns  pulse  shapes 
discussed  in  Section  5.3.2. 

5.4  Symmetric  removal  of  the  sidelobes  of  the  sine  function 

In  another  experiment,  we  symmetrically  remove  one  of  the  sidelobes  of  the  siiK  function, 
and  examine  the  effect  of  this  frequency  content  removal  on  the  shape  of  the  pulse.  We  use  a 
100  ns  pulse,  as  Figure  5. 15  shows,  for  these  experiments  so  that  the  sidelobes  are  10  MHz  wide, 
and  the  mainlobe  is  20  MHz  wide.  After  we  have  symmetrically  removed  the  sidelobes,  we 
perform  another  experiment  where  die  mainlobe  is  removed.  The  final  experiment,  presented  in 
this  section,  we  remove  all  the  sidelobes,  with  the  mainlobe  remaining  intact,  and  examine  the 
pulse  shape. 

Figure  5.16  shows  the  pulse  after  the  second  sidelobes  are  removed.  Notice  that  there  is  a 
sinusoid  on  the  passband  of  the  pulse  with  two  cycles  across  it  We  therefore  conclude  that  the 
second  sidelobes  provide  significant  energy  to  help  flatten  the  pulse.  This  pattern  continues,  as 
successive  sidelobes  are  removed,  to  generate  a  sinusoid  with  a  number  of  cycles  equal  to  the 
sidelobc  number  is  present  on  the  top  of  the  pulse.  For  example.  Figure  5.17  shows  the  pulse 
with  the  sixth  sidelobes  removed;  six  cycles  appear  across  the  top  of  the  pulse. 

Figure  5.18a  shows  that  when  we  remove  the  mainlobe  of  the  sine  function,  the  pulse  is 
devoid  of  the  d.c.  value  and  of  the  fundamental  The  remaining  sidelobe  spectral  energy  contains 
the  hi^er  frequency  components,  which  produce  the  sharp  rise-time  and  fall-time.  Figure  5.18b 
shows  the  pulse  with  all  the  sidelobes  of  the  sine  function  removed,  but  with  die  mainlobe  intact 
We  note  the  presence  of  the  d.c.  value  plus  the  fundamental  frequency,  which  contribute 
significantly  to  the  pulse  shape.  Figure  5.19a  and  5.19b  show  the  computer  simulation  results 
when  the  mainlobe  is  removed  and  then  all  the  sidelobes  of  the  sine  function  are  removed, 
respectively.  Notice  the  similarity  in  the  pulse  shapes  when  they  are  compared  to  Hgure  5. 1 8a  and 
5.18b.  Therefore  having  the  experimental  results  agree,  in  general,  with  the  computer  simulations. 


5.5  How  many  sidelobes  make  a  good  pulse? 

In  this  section  we  answer  the  question  “How  many  sidelobes  of  the  sine  function  in  the 
frequency  planers  needed  to  represent  a  good  pulse  in  the  time  domain?”.  In  this  e}q)enment  we 
use  a  100  ns  pulse  and  begin  the  experiment  with  the  spatial  li^t  modulator  elcmoits  switched  so 
that  seven  symmetric  sidelobes  are  present  Successive  sidelobes  are  truncated,  starting  at  the 
seventh  sidelobe,  and  proceeding  inward  towards  the  center  of  the  sine  function.  Figure  5.20a  and 
5.20b  show  the  computer  simulation  and  experimentally  obtained  pulse  when  seven  sidelobes  are 
present  in  the  frequency  plane.  From  the  data,  we  conclude  that  seven  sidelobes  clearly 
reconstructs  the  pulse. 

Figure  5.21  shows  the  experimental  pulse  when  five  sidelobes  are  present  in  the  frequency 
plane;  note  that  there  is  an  inaease  in  the  ripple  from  7  sidelobes  to  5  sidelobes.  The  rise-time  and 
fall-time  of  the  pulse  is  starting  to  increase  due  to  die  loss  of  the  frequency  conqionents  in  the  six 
and  seventh  sidelobes.  We  conclude  that  seven  tidelobes  is  the  minimum  to  cleanly  reconstruct  a 
pulse  and,  for  some  applications,  six  sidelobes  can  be  used. 

5.6  Introduction  to  narrowband  interference  experiments 

Figure  5.22  shows  the  setup  used  to  perform  narrowband  interference  excision 

experiments.  After  the  pulse  is  translated  to  f^,  a  narrowband  interference  signal  is  added,  the 

composite  signal  is  passed  through  the  Presort  Processor  and  is  demodulated  after  detectiotL  We 

denote  the  narrowband  interferer  as  j(t)  and  represent  the  signal  fed  into  the  Presort  Processor  as 

f(t)  =  s(t)+j(t).  (5.3) 

For  these  eiqieriments  the  narrowband  interference  is  a  single  frequency  at  the  frequency  fj. 

Therefore,  equation  (5.3)  becomes 

f(t)  *  s(t)cos(2n  fgt)  +  AjCOs(2it  fjt),  (5.4) 

where  f(t)  is  the  signal  into  the  acousto-optic  cell,  the  signal  is  Fourier  transformed,  the  appropriate 

spatial  light  modulator  elements  are  switched,  and  the  output  signal  is  translated  to  baseband  to 

recover  the  signal  s(t). 
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5.6.1  Excision  of  a  narrowband  interferer 

In  this  section  we  present  the  results  of  adding  a  narrowband  interferer  to  the  signal,  and 
perfonning  die  excision  of  the  interferer  to  illustrate  die  operation  of  the  Presort  Processor  The 
power  of  die  local  oscillator  is  adjusted  so  that  the  acousto-optic  cell  is  not  overdriven.  The 
reference  pulse  for  the  remainder  of  die  experiments  is  shown  in  Figure  5.23;  note  that  the  pulse 
has  a  5  mV  peak  and  a  ividth  of  100  nanoseconds.  Recall  that  the  sidelobes  are  10  MHz  wide  dius 
giving  a  mainlobe  width  of  20  MHz. 

A  narrowband  interferer  is  introduced  at  307  MHz,  which  is  in  the  coiter  of  the  second 
sidelobe  of  the  pulse.  Figure  524  shows  the  spectrum  of  the  narrowband  interferer  plus  the  pulse 
at  die  intermediate  frequency.  Notice  the  narrowband  interferer  is  30  dB  above  the  noise  floor. 
Figure  5.25  shows  the  narrowband  interferer  plus  the  pulse  waveform,  showing  that  the  pulse  is 
completely  obscured. 

Figure  5.26  shows  the  pulse  when  spatial  li^t  modulator  elements  centered  at  307  are 
switched  to  excise  the  narrowband  interferer.  The  output  pulse  has  the  same  shape  as  the  reference 
pulse  except  for  the  rise-time  overshoot,  which  is  apparent  in  the  reference  pulse,  that  is  partially 
removed  when  we  excise  the  narrowband  interferer.  This  result  illustrates  that  the  Presort 
Processor  is  capable  of  removing  narrowband  interferers,  without  seriously  afrecting  the  pulse 
shape. 

5.6.2  Notching  capability  of  the  Presort  Processor 

In  this  experiment  we  explore  the  notching  capability  of  the  Presort  Processor  by  locating 
the  narrowband  interferer  between  two  spatial  light  modulator  elements.  We  use  the  same  strength 
narrowband  interferer  as  in  the  first  experiment,  but  we  locate  the  narrowband  interferer  at  304.5 
MHz.  We  then  notch  only  two  spatial  light  modulator  elements  at  304  and  305,  with  Figure  5.27 
showing  the  output  pulse  shape.  Note  that  the  excision  of  the  narrowband  interferer  is  nearly 
conqilete,  and  a  comparison  with  the  reference  pulse  from  Fgure  5.23  reveals  a  good  shape  match. 
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The  differences  are  die  distortion  caused  by  excising  part  of  the  spectrum  of  the  pulse.  We 
dierefore  conclude  diat  if  a  narrowband  inteiferer  is  centered  between  two  spatial  li^t  modulator 
elements,  only  those  two  elements  need  be  switched 

5.7  Root  mean  square  error  of  pulse  after  a  narrowband  interferer  is  excised 

Since  90  percent  of  die  energy  is  in  the  mainlobe  of  a  sine  function  we  are  hi^y  concerned 
with  the  resulting  pulse  shape  when  the  narrowband  interferer  is  inside  the  mainlobe  of  die  sine 
function.  We  use  a  100  ns  pulse  which  gives  a  20  MHz  wide  mainlobe,  centered  at  331.7  MHz. 

Figure  5.28  shows  the  spectrum  of  the  pulse  plus  the  narrowband  interferec  Note  diat  die 
narrowband  interferer  is  roughly  30  dB  above  the  noise  level.  Suppose  that  we  move  the 
narrowband  interferer  through  the  mainlobe  from  left  to  right,  distorting  the  spectrum  of  the  pulse 
as  the  interferer  is  excised 

Figure  5.29  shows  the  root  mean  square  error  as  a  function  of  the  frequency  of  the 
interferec  The  root  mean  square  error  is  smallest  at  the  edge  of  the  mainlobe,  which  is  intuitively 
CQirea  since  the  mainlobe  has  a  lower  aiiqilitude  at  die  notch  position. 

As  we  increment  the  frequency  of  the  narrowband  interferer  through  the  mainlobe  we 
notice  that  the  root  mean  square  error  increases  until  332  MHz,  which  is  approximately  the  center 
frequency  of  the  pulse,  and  then  decreases  until  the  first  null  at  342  MHz.  Since  the  excision  of  the 
pulses  was  thoroughly  explored  in  earlier  sections,  we  refer  the  reader  to  those  sections  to  get  a 
flavor  of  the  pulse  shapes  with  these  frequency  components  excised 

As  a  final  example,  consider  the  case  where  the  narrowband  interferer  is  at  330  MHz. 

Since  this  frequency  is  close  to  the  center  frequency  of  the  pulse,  we  should  see  a  large  anplitude 
2  MHz  beat  frequency,  if  the  interferer  is  not  removed  Figure  5.30  shows  the  narrowband 
interferer  obscuring  the  pulse  waveform.  Figure  5.31  shows  the  pulse  when  we  form  the 
appropriate  notch  at  330.  This  is  an  extremely  pleasing  and  dynamic  result  The  pulse  waveform 
is  largely  preserved,  with  some  distortion  caused  by  the  removal  of  some  frequency  components. 


6.  Conclusions 
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Frequency  plane  excision  is  one  method  for  removing  narrowband  interference  from 
wideband  signals.  We  showed  that  the  Presort  Processor,  which  utilizes  frequency  plane  excision, 
can  be  modified  so  that  the  fiequency  excision  can  evenmally  be  automated.  The  control  signals  to 
the  spatial  light  modulators,  which  perform  die  notching  of  the  interference,  are  updated  every  5.4 
milliseconds  exceeding  the  design  goal  of  10  milliseconds.  We  showed  both  experimentally  and 
dieoretically  that  the  proposed  modification  of  the  Preson  Processor  did  not  affect  its  performance. 

We  analyzed  area  detection  in  the  Fourier  plane  versus  point  detection  in  die  image  plane 
and  showed  that  the  two  detection  schemes  are  essentially  equal  and  produce  notches  of  equal 
depth.  We  further  showed  that  the  point  detection  signal  bias  tenn  is  slightly  less  than  the  area 
detection  signal  bias  term,  thereby  reducing  the  shot-noise.  We  also  developed  an  image  plane 
analysis  technique  that  adds  insight  into  the  into'-reladonship  of  the  parameters  used  in  fiequency 
plane  excision. 

We  tested  the  operation  of  the  Presort  Processor  by  removing  various  fiequency 
components  of  a  pulse  and  analyzing  the  resulting  distortion  at  baseband.  As  expected,  these 
experiments  showed  that  the  distortion  is  most  severe  when  the  removal  is  performed  in  die 
mainlobe  of  the  sine  function.  Further  experiments  added  insight  into  the  relationship  between 
removing  symmetric  sidelobes  in  the  fiequency  plane  and  the  resulting  distortion  of  the  image 
plane  pulse.  We  also  added  a  narrowband  inierferer  to  the  pulse,  performed  the  excision  using  the 
Presort  Processor,  and  analyzed  the  image  plane  pulse.  These  experiments  showed,  as  expected, 
that  the  distortion  on  the  pulse  is  small  if  the  interferer  is  outside  the  mainlobe  and  is  more 
pronounced  if  the  interferer  is  inside  the  mainlobe  of  the  sine  function. 

Suggestions  for  future  research  include: 

•  increasing  the  speed  of  the  post-processing  electronics  by  using  higher  speed  logic  chips. 

•  performing  multi-interference  experiments  where  the  interferers  can  be  agile  or  stationary. 
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Figure  2.2:  Optical  Fouiio’  transform  system 
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Figure  2.3:  Parameters  of  an  acousto-optic  cell 
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Figure  2.6;  Relationships  of  the  acousto-optic  c 
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Figure  2,7:  Fourier  transform  configuration 


amplitude 


ampUtude 


Figure  4,1;  Preson  Processor 


Figure  4.2:  Passband  of  the  signal  acousto-opdc  cell 
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Figiire  4.7:  Amplitude  and  phase  response  of  the  Preson  Processor 


(a)  phase  response  after  window  is  insened 


(b)  phase  response  after  linear  phase  is  nulled 
Figure  4.8:  Phase  response  after  the  optical  window  is  inserted  into  the  Presort  Processor 
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(a)  Expanded  view  of  the  amplitude  response  before  the  optical  window  is  inserted 
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(b)  Expanded  view  of  the  amplitude  response  after  the  optical  window  is  inserted 
Figure  4.9:  Amplitude  response  before  and  after  the  optical  window  is  inserted 
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Figure  4. 10;  Phase  response  after  the  corrective  optical  window  is  positioned  in  the  upper  branch 


Figure  4.11;  Pulse  before  the  primary  optical  window  is  inserted,  after  the  pnmary  optical 
window  is  inserted,  and  after  the  corrective  optical  window  is  inserted. 
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Figure  4.14:  Block  diagram  of  the  digital  system 
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Figure  4.15:  Specmun  of  a  signal  distorted  by  narrowband  interferers 


84 


Rgurr  4. 16;  The  display  with  one  nanowband  interferer  excised 
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Bgure  4.17;  Experimental  setup  for  the  testing  of  the  post-processing  electronics 
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Hgure  4. 18:  The  spectrum  with  the  threshold  level  at  a  steady- state  value 


Figure  4.19:  The  spectrum  with  with  two  different  strength  interfercrs  with  the  threshold  at  a 
steady  state  value 


86 


Figure  4.20:  The  spectrum  with  multiple  narrowband  interferers 
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Figure  5.1:  Block  diagram  of  the  experimental  system 


Hgure  5.2:  Experimental  equipment  setup 


Figure  5.3:  Signal  that  is  available  on  a  one-shot  basis 


Figure  5.5:  Computer  simulation:  root  mean  square  error  versus  spatial  light  modulators  elements 


Hguie  5.6:  Experimental  results:  root  mean  square  error  versus  spatial  light  modulators  elements 
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Figure  5.9:  The  78  ns  pulse  with  the  d.c.  value  removed 


(a)  Computer  simulation  result:  spatial  light  modulator  elements  345  and  355  switched 


(a)  Computer  simulation  result:  spatial  light  modulator  elements  332  and  368  switched 
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(a)  Computer  simulation:  root  mean  square  error  versus  spatial  light  modulators  elements 
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(b)  Experimental  results:  root  mean  square  error  versus  spatial  light  modulators  elements 
Figure  5.13:  Root  mean  square  error  of  the  468  ns  pulse 


(a)  The  100  ns  pulse  with  the  mainlobe  removed 


(a)  Computer  simulation  result  of  the  100  ns  pulse  with  the  mainlobe  removed 
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(b)  Computer  simulation  result  of  the  100  ns  pulse  with  the  sidelobes  removed 
Figure  5.19:  The  computer  simulations  of  the  100  ns  pulse  with  various  sidelobes  removed 


(a)  Computer  simulation  result  of  the  100  ns  pulse  with  seven  sidelobes 
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Figure  5.22:  The  experimental  setup  for  the  narrowband  interference  experiments 
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Figure  5.23:  The  reference  pulse  for  the  nanowband  interference  experiments 
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Figijre  5.24:  The  spectrum  of  the  pulse  plus  the  narrowband  interfener 
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Figure  5.25:  The  pulse  plus  the  narrowband  interfercr  after  demodulation 
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Figure  5.30;  The  pulse  plus  the  narrowband  interferer  2  MHz  from  the  modulating  frequency 
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Figure  5.31:  The  pulse  with  the  narrowband  interferer  excised:  RMSE=0.15 
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ROME  LABORATORY 

Rome  Laboratory  plans  and  executes  an  interdisciplinary  program  in  re¬ 
search,  development,  test,  and  technology  transition  in  support  of  Air 
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Force  Command,  Control,  Communications  and  Intelligence  (C  I)  activities 
for  all  Air  Force  platforms.  It  also  executes  selected  acquisition  programs 
in  several  areas  of  expertise.  Technical  and  engineering  support  within 
areas  of  competence  is  provided  to  ESD  Program  Offices  (POs)  and  other 
ESD  elements  to  perform  effective  acquisition  of  C  I  systems.  In  addition, 
Rome  Laboratory's  technology  supports  other  AFSC  Product  Divisions,  the 
Air  Force  user  community,  and  other  DOD  and  non-DOD  agencies.  Rome 
Laboratory  maintains  technical  competence  and  research  programs  in  areas 
including,  but  not  limited  to,  communications,  command  and  control,  battle 
management,  intelligence  information  processing,  computational  sciences 
and  software  producibility,  wide  area  surveillance/sensors,  signal  proces¬ 
sing,  solid  state  sciences,  photonics,  electromagnetic  technology,  super¬ 
conductivity,  and  electronic  reliability/maintainability  and  testability. 


